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Abstract 

This  paper  models  the  absorption  phenomena  of  light  in  a  three  level  diode 
pumped  alkali  laser  system.  Specifically  this  model  calculates  for  a  user  defined  set  of 
system  parameters  the  attenuation  of  the  input  pump  beam  and  characteristics  of  the 
bleached  wave.  Using  Wolfram’s  Mathematica  6.0  software  all  necessary  physics  for  an 
accurate  description  of  absorption  was  modeled  from  first  principles:  energy  levels,  cross 
sections,  spin-orbit  kinetic  processes,  saturation  frequencies,  pump  attenuation,  and 
differential  transmittance,  which  is  a  representation  of  the  bleached  wave.  A  specific 
DPAL  scenario  was  simulated,  455K  system  temperature,  alkali  concentration  of  6.1  • 

10  ,  and  a  system  pressure  of  200  torr  of  He  and  600  torr  of  Ethane.  For  a  range  of  initial 
input  intensities  the  linear  approximation  to  the  beam  attenuation  predicted  consistently  a 
differential  transmittance  value  of  70%.  It  was  concluded  that  the  linear  approximation  is 
a  good  indicator  of  the  distance  a  bleached  wave  penetrates  an  absorbing  cell.  This  model 
was  also  benched  marked  against  the  quasi  two  level  model.  In  the  limit  of  high  system 
pressure  the  simulated  model  converged  to  the  same  population  inversion  as  that  of  the 
quasi  two  level  regime.  Finally,  within  the  quasi  two  level  regime  a  closed  analytic 
equation  was  developed  to  describe  under  what  conditions  the  system  would  lase. 
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A  THEORETICAL  MODEL  ANALYSIS  ON  ABSORPTION  OF  A  THREE 
LEVEL  DIODE  PUMPED  ALKALI  LASER 

I.  Introduction 

The  purpose  of  this  model  is  to  explore  the  key  aspects  of  absorption  in  a  diode 
pumped  alkali  lasing  (DPAL)  system.  The  key  performance  parameter  space  will  be 
explored  and  reduced  to  controllable  aspects  that  will  aid  in  the  engineering  of  such  a 
system.  To  date  the  United  States  Air  Force  has  invested  heavily  in  chemical  lasing 
systems.  Although  lasing  powers  have  been  achieved  to  become  both  tactically  and 
strategically  valuable  (megawatt  class),  the  logistic  footprint  and  finite  magazines  of 
deploying  chemical  lasing  systems  have  become  a  cause  to  research  other  lasing  systems. 
Research  has  since  been  turned  toward  achieving  power  levels  on  the  order  of  lOOkW 
with  solid  state  lasers.  Progress  in  attaining  such  power  levels  have  not  been  as  fast  as 
some  researchers  have  hoped.  As  a  possible  solution  to  this  problem  the  Air  Force 
Institute  of  Technology  along  with  its  technology  partners  have  been  tasked  to  see  if  a 
combination  laser,  one  that  is  part  solid  state  and  part  gas,  could  fit  the  requirements  of 
being  both  cost  effective  and  tactically  viable.  Many  aspects  must  come  together  in  order 
for  this  to  happen.  To  date  most  research  has  been  done  experimentally  through  trial  and 
error. 

A  DPAL  system  is  in  essence  a  gas  electric  hybrid  laser.  The  lasing  medium  is  a 
gaseous  cell  of  alkali  typically  Rubidium,  Cesium,  or  Sodium.  The  unique  character  of 
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alkali  atoms,  a  single  valence  electron,  allows  its  first  two  excited  states  to  be  taken 
advantage  of.  A  diode  electric  laser  is  used  to  pump  or  excite  the  alkali  atoms  to  its 
second  excited  state  and  then  inert  gases  are  introduced  in  to  the  cell  to  collisionally  de- 
excite  the  alkali  atoms  to  its  first  excited  state.  The  energy  released  then  as  the  atoms 
finally  return  to  their  ground  state  is  in  the  form  of  coherent  lasing  light.  One  big 
advantage  of  this  system  is  thermal  control.  The  flow  rate  of  the  gaseous  medium  through 
the  lasing  cell  can  be  controlled  thereby  at  the  same  time  removing  heat  from  the  DPAL 
system.  However  there  exists  the  technical  challenge  of  matching  the  bandwidth  of  the 
diode  laser  stacks  to  that  of  the  absorption  bandwidth  of  the  alkali  medium.  This  can  be 
slightly  alleviated  by  broadening  the  alkali  absorption  spectrum  by  running  the  cell  at  a 
high  pressure. 

It  is  the  aim  of  this  paper  to  create  a  theoretical  tool  to  tackle  one  aspect  of  a 
DPAL  laser,  which  is  the  absorption  of  light  in  the  gas  lasing  medium.  The  problem  is  to 
first  describe  the  absorption  cross  section  of  Cesium  and  Rubidium  to  the  hyperfine 
structure  level  under  various  system  temperatures  and  pressures.  Second,  is  to  solve  for 
the  rate  at  which  populations  in  the  three  different  energy  levels  change  under  various 
system  conditions.  This  will  allow  the  deduction  of  how  the  diode  pump  field  photons  are 
absorbed  or  attenuated  by  the  alkali  atoms  under  both  saturated  and  unsaturated 
conditions.  From  these  ideas  a  picture  of  what  a  bleached  wave  is  and  how  it  affects 
system  absorbance  will  be  explored.  Next,  these  ideas  will  be  extended  to  a  quasi-two 
level  model  in  which  the  two  excited  states  are  infinitely  coupled  or  the  system  is 
operating  at  very  high  pressure.  The  last  challenge  will  be  to  introduce  a  lasing  field  into 
the  quasi-two  level  regime  and  derive  an  analytic  equation  that  will  determine  if  under 
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what  system  conditions  is  lasing  possible.  This  model  will  be  simulated  in  a  computing 
environment  to  help  DPAL  researchers  determine  what  parameters  a  full  up  system  will 
most  likely  be  sensitive  to  as  well  as  determine  certain  configurations  that  would  limit 
absorption  or  lasing  performance.  This  is  not  intended  to  be  a  full  model  of  a  tactically 
capable  system,  but  rather  a  component  of  what  could  become  a  full  up  model. 
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II.  Theory 


Chapter  Overview 

The  theory  chapter  will  be  devoted  to  the  explanation  and  derivation  of  the 
necessary  physical  theories  that  this  model  will  use  to  simulate  the  absorption  of  light  by 
alkali  atoms.  The  theories  contained  in  this  chapter  are  chiefly  used  by  atomic 
spectroscopist  and  rely  heavily  on  physics  of  light  and  matter  interaction  as  well  as 
quantum  mechanics  and  kinetics.  First,  the  properties  and  structure  of  alkali  atoms  will  be 
discussed  to  the  necessary  precision  required  of  this  model.  Next,  how  light  interacts  with 
matter  will  be  discussed  to  understand  how  and  why  a  gas  of  alkali  atoms  would  absorb 
incoming  photons.  Then,  a  simple  discussion  of  how  the  systems  of  alkali  atoms  change 
behavior  when  foreign  species,  other  atoms  or  molecules,  interfere  with  the  absorption  of 
light.  Lastly,  to  mimic  real  conditions  the  effects  of  uncertainty,  temperature,  and 
pressure  must  be  dealt  with  if  this  model  is  going  to  be  useful. 

Alkali  Atomic  Structure 

Properties  and  Energy  Levels 

According  to  the  International  Union  of  Pure  and  Applied  Chemistry  (IUPAC)  the 
series  of  alkali  metal  atoms,  lithium,  sodium,  potassium,  rubidium,  cesium,  and  francium, 
comprise  of  the  first  family  of  elements  on  the  periodic  table.  What  makes  these  elements 
chemically  similar  and  which  defines  their  family  is  the  position  of  their  valence  electron: 
it  solely  occupies  an  S-orbital.  This  observation  is  critical  to  how  the  alkali  metals 
interact  with  their  surrounding  environment.  Table  1  below  lists  some  basic  physical 
properties  of  the  alkali  metals: 
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Table  1.  Alkali  Physical  Properties 


Lithium 

Sodium 

Potassium 

Rubidium 

Cesium 

Francium 

Atomic  Number 

3 

n 

19 

37 

55 

87 

Atomic  Weight 
(amu) 

6.941 

22.990 

39.098 

85.468 

132.905 

223 

Melting  Point 

(°C) 

180.54 

97.72 

63.38 

39.31 

28.44 

Unknown 

Boiling  Point 
(°C) 

1342 

883 

759 

688 

671 

Unknown 

Specific  Heat 
Capacity 
(J/mol-K) 

24.860 

28.230 

29.597 

31.06 

32.210 

Unknown 

Nuclear  Spin 

3/2 

3/2 

3/2 

3/2(87) 

7/2  (85) 

7/2 

3/2 

Ionization  Limit 
(kJ/mol) 

520.2 

495.8 

418.8 

403.0 

375.7 

Unknown 

The  similarity  of  these  alkali  atoms  to  the  simplest  atom,  hydrogen,  makes  it  easy 
to  understand  the  energy  levels  of  these  atoms.  The  hydrogen  atom  consists  of  one 
electron  orbiting  one  proton.  Likewise  the  alkali  metals  consist  of  a  positively  charged 
core,  the  protons  in  the  nucleus  plus  the  inner  electrons  add  to  +1  charge,  with  a  single 
electron  inhabiting  the  outermost  orbit.  First,  start  with  the  nonrelativistic  Hamiltonian 
operator  describing  the  energy  of  the  valence  electron  for  the  Hydrogen  Atom: 


H  = 


(1.1) 


where, 

h  =  Planck’s  constant 

me  =  mass  of  the  electron  (kg) 

e  =  charge  of  the  electron  (coulomb) 

r  =  distance  separating  proton  and  electron  (m) 
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The  electron  experiences  the  familiar  Coulomb  potential  which  in  this  case  is  also 
spherically  symmetrical.  Second,  cast  the  operator  into  spherical-polar  coordinates.  The 
time  independent  Schrodinger  equation,  HXP=E'¥ ,  can  now  be  solved  using  the  common 
separation  of  variables  method.  The  solution  for  both  the  energy  eigenstates  and  the 
normalized  wavefunctions  are  as  follows: 


~mee 
lh2n 2 


(1.2) 


a\n-l- 1)! 
2«[(«  +  l)!] 


_p. 


e-'VC'(p)r(M) 


(1.3) 


where, 


n  =  electronic  energy  level  quantum  number 
1,  m  =  rotational  energy  level  quantum  numbers 
2 

a  = - r ,  a0  =  Bohr  Radius 

na0 

2  r 

P  = - 

na0 

C/(P)  :  Associated  Laguerre  Polynomials 
Y"‘(0,(j))  :  Normalized  Spherical  Harmonics 
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The  probability  distributions  of  these  wavefunctions  give  rise  to  the  familiar  electron 
orbitals  learned  in  basic  chemistry,  s,  p,  d,  f,  etc. 

The  alkali  metal  elements  can  be  cast  into  the  same  analytical  ansatz  as  the 
hydrogen  atom.  The  outermost  electron  of  the  alkalis  experience  the  same  spherical 
Coulomb  potential  as  that  of  hydrogen.  As  stated  before  adding  both  the  protons  and 
electrons,  excluding  the  valence  electron,  add  up  to  a  charge  of  +1.  Noticing  this  trend  the 
energy  levels  follow  the  same  functional  form,  eqn.  1.2,  as  that  of  hydrogen  with  the 
following  modifications: 


EnJ  = 


-m/Z2 
2h2  (n  —  8n  jj2 


(1.4) 


where,  Z  =  Z-  N+  l,Zis  the  number  of  protons  in  the  nucleus,  and  N  is  the  number  of 
electrons.  The  parameter  8n ,  contains  corrections  factors  that  account  for  the  shift  from 

hydrogen  behavior.  The  first  correction  is  a(l) ,  the  quantum  defect  factor,  which  is  a 
function  of  the  valence  electron’s  orbital  momentum  value.  As  the  electron  orbits  the  core 
of  protons  and  electrons  it  is  quite  possible  that  as  higher  angular  momentum  states  are 
achieved  the  distance  in  which  the  electron  penetrates  the  core  will  be  greater.  This  leads 
to  the  second  correction  factor  f3(l)  which  accounts  for  the  fact  that  as  the  electron 
burrows  into  the  core  it  experiences  less  shielding  from  the  core  electrons  and  so  the 
Coulomb  potential  looks  less  classical.  This  is  particularly  important  the  closer  the 
valence  electron  gets  to  the  nucleus.  The  final  form  of  this  correction  parameter  is, 
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8nj=a(l)  + 


(1.5) 


1 3(1) 


n 


2 


Table  2.  Quantum  defect  parameters  for  Li,  Na,  K,  Rb,  and  Cs 


1  = 

0 

1 

2 

3 

Li 

0.40 

0.04 

0.00 

0.00 

Na 

1.35 

0.85 

0.01 

0.00 

K 

2.19 

1.71 

0.25 

0.00 

Rb 

3.13 

2.66 

1.34 

0.01 

Cs 

4.06 

3.59 

2.46 

0.02 

Fine  Structure  Splitting 

Another  important  phenomenon  for  the  complete  picture  of  the  alkali  atom 
structure  is  the  coupling  of  both  the  spin  and  orbital  angular  momentum  of  the  electron. 
For  the  hydrogen  atom  this  spin-orbit  (SO)  coupling  can  be  understood  in  a  semi- 
classically  way  by  noting  that  in  the  electron’s  rest  frame  the  proton  is  orbiting  the 
electron  and  this  movement  leads  to  current.  This  current  in  turn  invokes  a  magnetic  field 
and  the  electron’s  spin  then  couples  to  this  field.  However,  since  SO  coupling  is 
inherently  a  relativistic  phenomenon  it  is  natural  to  start  with  the  Dirac  equation  that  is 
appropriate  for  spin  “A  particles  i.e.  electrons. 

Start  with  the  classical  relativistic  energy  as  prescribed  by  Einstein: 


=  {c2  p2 


+  m0  c 


\l/2 


(1.6) 


where, 


c  =  speed  of  light  (m/sec) 


8 


p  =  momentum  (kg-m/sec) 
m0  =  rest  mass  (kg) 

When  this  Hamiltonian  is  substituted  into  the  Schrodinger  equation  and  the  square  root 
expanded  in  a  series  of  the  momentum,  it  treats  time  and  space  asymmetrically  which 
cannot  be  the  case  since  time  and  space  are  symmetric  in  special  relativity.  Therefore  a 
modification  must  be  performed  first  so  that  they  become  symmetric  again. 

This  modification  involves  transforming  the  term  in  the  parentheses  of  equation 
1.6  to  a  perfect  square.  In  this  way  the  Schrodinger  equation  will  be  both  linear  in  space 
and  time: 


c2P2  +  m2cA  =  (. caxPx  +  cavP  +  cazP  +  pm0c2  f 
=  (ca-P  +  / 3m0c 2)2 

The  constraints  imposed  by  this  method  give  the  following  requirements  for  a  and  P: 

a2=p2=  1  (i  =  x,y,z) 

aiaj  +  ajai  -  [a, ,  aj  ]  =  0  (/'  ^  j)  (1.8) 

aip  +  pai  =  [ai,p\=  0 


In  order  for  H  to  remain  Hermitian  the  a,  p  matrices  must  be  themselves  Hermitian  as 
well  as  being  traceless  and  having  eigenvalues  of  ±1.  These  three  requirements  mean  that 
the  matrices  that  complete  the  perfect  square  of  the  relativistic  Hamiltonian  must  be  4x4. 
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a  = 


0  o 
o  0 


P  = 


I  0 

0  -I 


(1.9) 


The  choice  of  what  a,  [S  are  not  unique  since  they  can  be  transformed  by  a  matrix  that  is 
unitary.  The  a  are  the  familiar  2x2  Pauli  matrices  and  the  identity  matrices  are  2x2.  With 
such  parameters  defined  Paul  Dirac’s  equation  can  be  fully  realized: 


iti 


jM 

dt 


(ca-P  +  pw0c2)|t//^ 


(1.10) 


The  state  vectors  are  now  a  four-component  object  rather  than  the  two-component  object 
familiar  in  non-relativistic  quantum  mechanics.  However  the  energy  eigenstates  look 
vastly  familiar  with  the  modification  that  it  contains  two  two-compnent  spinors,  %  and  <f>: 


V(0  = 


X 

O 


g-iEt/h 


(1.11) 


Using  the  well-known  time-independent  Schrodinger  equation  the  spin  orbit 
interaction  will  naturally  fall  out  as  long  as  orders  of  the  equation  are  expanded  out  to 

(v/  c)4 .  First,  start  by  inserting  the  Coulomb  potential,  V  =  —e1  /  r ,  into  the  relativistic 

Hamiltonian  which  will  yield  a  set  of  coupled  spinor  equations.  The  first  spinor  will 
describe  the  simple  case  of  an  electron  orbiting  a  fixed  proton.  The  second  spinor 
contains  information  about  a  positron  which  is  not  of  interest  in  this  paper. 
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(1.12) 


E-V  -  m0c2 

—CG  •  P 

~X~ 

-ca-  P 

E-V  -  m0c2 

o 

Solving  for  O  and  then  inserting  into  the  equation  for  %  will  yield  an  uncoupled  equation 
in  x  only: 


(E  -V  -  m0c2  )x  =  cg-  P 


1 


E-V  +  mnc~ 


ca-Px 


(1.13) 


Now  a  couple  of  approximations  involving  the  total  energy,  both  kinetic  and  rest,  can  be 
made  to  recover  the  non-relativistic  Schrodinger  equation.  Define  the  Schrodinger 
energy,  Es,  as  the  total  energy  minus  the  rest  energy.  If  we  add  one  unit  of  rest  mass 
energy  to  the  total  energy  this  relation  appears: 


E  +  in qC  -  Es+  2 m0c 


(1.14) 


At  very  low  speeds  compared  to  the  speed  of  light  the  Schrodinger  energy  is  much  less 
than  the  rest  energy  and  so: 


E  +  m0c2  =  2m0c2 


(1.15) 
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2 

This  allows  the  denominator  in  equation  1.13  to  be  approximated  as  2moc~  assuming  the 
potential  doesn’t  make  any  appreciable  contributions.  Solving  for  the  electron  spinor  the 
standard  Schrodinger  equation  is  recovered. 


E,X  = 


(?2l 

2  mn 


+  V 


X 


2  mn 


-  +  V 


X 


(1.16) 


However  in  order  to  see  the  fine  structure  the  denominator  in  equation  1.13  has  to  be 
expanded  to  P4  or  (v/c)4. 


1  1  1  1 

fi  Es-V) 

I  I  S 

E-V  +  m0c2  Es-V 

1 

+  2  m0c2  2 

m0c 2 

IT  9 

^  2m0c~  J 

1  E-V 

2  m0c2 

2  mnc2 

\  0  J 

2  m0c2  4  mac4 

(1.17) 


Inserting  this  back  into  equation  1.13  it  begins  to  look  closer  to  an  analyzable 
Hamiltonian: 


ESX  = 


p2  ¥V  o-P(£,-r)g-P 
2  m0  4  mgC2 


X 


(1.18) 
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If  the  quantity  Es  -  V  acted  directly  on  y  the  first  order  solution  could  simply  be 
substituted.  However  since  this  is  not  the  case  a  few  manipulations  using  the  commutator 
can  rectify  the  situation. 


(es-v)*y>x=«v(es-v)x+g\es-v,v]x 

P2  (1-19) 

=  0-p)— x+[v,v]x 

2m0 


Equation  1.13  now  becomes  and  equation  in  which  the  Hamiltonian  has  two  familiar 
terms  and  three  terms  due  solely  to  relativity: 


E,X  = 


2  mn 


-  +  V- 


/o  Px[P/]  P-[P,k]l 


8  m03c2 


4  m02c2 


4  m02c2 


(1.20) 


The  third  term  is  the  relativistic  correction  to  the  kinetic  energy,  the  fourth  term  contains 
the  correction  that  is  sought  after,  the  SO  coupling  energy,  and  the  fifth  term  is  known  as 
the  Darwin  correction.  The  evaluation  of  the  SO  coupling  term  is  straight  forward: 
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For  the  Alkali  metals  this  interaction  causes  a  splitting  of  the  energy  levels.  To 
compute  the  energy  corrections  one  has  to  move  to  a  basis  that  will  completely 
diagonalize  H.  In  this  case  the  total  angular  momentum  basis,  J=L+S,  is  sufficient  to 
understand  what  happens  at  the  fine  structure  level.  The  energy  corrections  are 
proportional  to  the  quantity: 


Eso  -  ^)-[./(/  +  l)-Z(T  +  l)-lS0S  +  l)]  (1.22) 


The  valence  electrons  in  each  Alkali  atom  in  the  ground  state  are  contained  in  an  s- 
orbital.  So  their  angular  momentum  numbers  are  J  =  1/2,  L  =  0,  and  S  =  ±1/2.  As  such  no 
fine  splitting  occurs  in  the  ground  state.  Once  the  electron  is  promoted  to  the  next  energy 
level  it  finds  itself  in  a  p-orbital  with  momentum  numbers  J  =  3/2,  1/2,  L  =  1,  and  S  = 
±1/2.  So  an  Alkali  valence  electron  when  excited  to  the  next  energy  level  has  two  fine 
structure  levels  to  ascend  towards.  This  gives  rise  to  the  characteristic  doublet  absorption 
lines  seen  in  the  alkali  metals. 
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Hyperfine  Structure  Splitting 

The  hyperfine  structure  will  be  the  furthest  this  paper  needs  to  go  into  the  internal 
structure  of  the  atom.  Qualitatively  the  hyperfine  structure  is  brought  about  by  the  same 
scheme  the  SO  correction  came  out  of:  interaction  of  electromagnetic  fields  and  particle 
spin.  The  fine  structure  coupled  the  spin  of  the  electron  with  its  orbital  momentum  via  an 
electromagnetic  interaction.  However  just  as  the  electron  has  intrinsic  spin  so  do  the 
particles  that  makes  up  the  nucleus.  This  is  called  nuclear  spin  and  is  typically  labeled  I. 
The  nucleus  can  be  attributed  a  magnetic  moment  and  therefore  depending  on  its  size  will 
have  multi-pole  moments  describing  its  field.  So  the  picture  looks  as  such  for  the 
hyperfine  interaction: 

1. )  The  magnetic  moment  of  the  nucleus  will  interact  with  the  magnetic  field 

produced  by  the  current  of  the  orbiting  electron. 

2. )  Both  the  nucleus  and  electron  by  virtue  of  their  spin  produce  magnetic 

moments.  As  such  multi-pole  interactions  occur. 

3. )  Lastly,  there  is  what  is  known  as  Fermi’s  contact  term  which  describes  the 

interaction  of  the  electrons  magnetic  moment  with  the  magnetic  field  inside 

the  nucleus  to  include  singularities  that  exist  inside  the  particles. 

Once  again  we  will  take  the  simple  case  of  the  hydrogen  atom,  determine  the  hyperfine 
Hamiltonian  from  first  principles,  and  then  determine  what  it  will  take  to  extend  to  alkali 
systems. 

The  Hamiltonian  of  an  electron  in  the  presence  of  a  magnetic  field  produced  by  a 
proton  is: 
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(1.23) 


H 


=  2m~ ~~ ^ (f )]2  ~~ SVl (? ) ~~  2Mb 


VxAj(r) 


where, 


Ai  =  Vector  potential  of  nucleus 

Pb  =  The  Bohr  Magneton 

Vi  =  electric  potential  of  the  nucleus 

Note  that  the  final  term  is  the  coupling  of  the  electron’s  spin  with  the  magnetic  field  of 
the  nucleus.  The  vector  potential  is  related  to  the  magnetic  dipole  moment  by: 


where, 


Mr) 


Mo  MTxf 
An  r 3 


Mi  =  magnetic  dipole  moment  of  the  nucleus 
p0  =  magnetic  permeability 


(1.24) 


Inserting  equation  1.24  into  equation  1.23  and  ignoring  any  quadratic  terms  of  the  vector 
potential  the  first  term  in  1.23  becomes: 


P-AI(r)  +  A,(f)-P 


Mo  [P-(Mtxf)  (Mtxr)-P 

An  r3  r3 


(1.25) 


The  magnetic  moment  commutes  with  both  P  and  R  since  it  only  involves  the  spin 
operator  so  the  first  and  second  terms  of  1.25  can  be  rewritten  as: 
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(1.26) 


P  (MjXr)_  L  Mj 

3  _  3 

r  r 

(M1xr)-P_M,L 

3  _  3 

r  r 

One  just  has  to  be  careful  when  using  cross  product  identities  not  to  upset  the  order  of 
multiplication  between  P  and  R  since  they  do  not  commute  like  the  scalar  operator  Irl.  As 
such  the  first  two  terms  in  equation  1.23  become  something  that  is  analyzable  by 
inspection: 


Bo  g  2  ^ 

4k  2me  r3 


-2M,  Bl 


(1.27) 


Now  it  is  clear  that  the  magnetic  moment  of  the  nucleus  is  coupling  to  the  magnetic  field 
due  to  the  current  of  the  electron. 

The  last  two  interactions,  the  dipole-dipole  and  Fermi  interaction,  are  contained  in 
the  last  term  of  equation  1.23.  This  involves  knowing  the  magnetic  field  both  outside  and 
inside  the  proton.  For  the  field  outside  the  field  is  simply  that  of  a  dipole  field.  If  we  align 
the  magnetic  dipole  moment  in  the  z-direction  the  dipole  field  is: 


x  4  n  r 5 


fl 


:^3 
4k  r 


B  = 


Bo_ 

4k 


M, 


3 z2  -  r 2 


(1.28) 
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Now  to  complete  the  dipole  interaction  term  we  take  the  components  of  the  electronic 
spin  operator  and  dot  product  them  into  equation  1.28.  Also  by  definition  Mi  is  in  the  z- 
direction  so: 


Wdip  - 


4  n 


¥ 


■Vxdy(r)  = 


V  J 


Mo  ^Mb_ 
4  /r/zr3 


S-Mj-3 


(Sf)(M,.S) 


(1.29) 


The  electronic  spin  is  now  coupled  to  the  nuclear  spin,  in  this  case  the  proton. 

The  magnetic  field  inside  the  proton  is  calculated  by  adding  the  flux  outside  the 
radius  of  the  proton,  po,  the  flux  inside  the  radius,  setting  them  equal  to  zero  and  solving 
for  B;: 


4k  p0 


(1.30) 


To  understand  the  effects  of  this  term  insert  1.30  into  the  last  term  of  equation  1.23  and 
use  the  |  n,l,m,s^  representation  to  expand: 


n,l, m ,s |  MiAMi’M i  | n\l\m\s'Y 


4ntip0 


Mo  ^Mb^ i 
Anhp* 


sz  |  *  ')  JJJ£p  ^r(p*n  l  m  (r)<p„, m,(r) 


(1.31) 
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4  , 

As  po  goes  to  zero  so  does  the  volume  element,  —  7rp0  >  and  so  the  integral  becomes  a 
delta  function  at  r  =  0.  Thus  the  last  hyperfine  correction: 


H 


Fermi 

¥ 


~Mo  Sx 
4  n  3 


M, 


'ln.S'' 


8(r) 


(1.32) 


This  completes  what  is  known  as  the  hyperfine  structure  of  the  Hydrogen  atom: 


H_  1 1 '  ,  tt  dip  .  -fj  Fermi 

¥  ~  **¥  +  ** /'/  +  **¥' 


(1.33) 


What  makes  this  analysis  so  straight  forward  is  that  the  nucleus  contains  only  one  proton 
and  so  only  the  dipole  term  needs  to  be  considered.  The  number  of  multipoles  that  must 
be  included  in  the  above  analysis  must  be  less  than  or  equal  to  twice  the  nuclear  spin. 
Equation  1.23  would  then  contain  extra  terms  in  both  the  electric  potential  and  magnetic 
potential.  For  the  alkali  metals  this  is  definitely  the  case  and  though  the  derivation  will 
not  be  repeated  the  same  process  will  yield  the  hyperfine  Hamiltonian  for  alkali  metals 
(using  a  suitable  angular  momentum  basis): 


H/rf-j  —  A  hJ  •  J  +  B 


3(1  •  J)2  +  2  (I  ■  J) -/(/  +  !)/( J  +  l) 


+ 


C 


hfs  hff  hfs  2/(2/  -  Y)J(2J  - 1) 

10(1  •  J)3  +  20(1  •  J)2  +  2(1  •  J)  [1(1  + 1)  +  J(J  + 1)  +  3]-  31(1  + 1)  J(J  + 1)  -  51(1  + 1  )J(J  + 1) 


hf  s 


I  (I  - 1)  (21  - 1  )J(J  - 1)  (2J  - 1) 


(1.34) 
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where  Ahfs  is  the  magnetic  dipole  constant,  Bhfs  is  the  electric  quadrupole  constant,  and 
Chfs  is  the  magnetic  octupole  moment.  Below  in  Table  x  are  the  multipole  moments  of 
Cesium  and  Rubidium. 


Table  3.  Hyperfine  Multipole  constants  for  Cesium  and  Rubidium 


Moment 

Cesium 

Rubidium  85 

Rubidium  87 

A  Sm 

h-2.2981  GHz 

h  l.0119  GHz 

h-3.4173  GHz 

AP1/2 

h-291.9201  MHz 

h- 120.527(6)  MHz 

h-408. 328(2)  MHz 

A  P3/2 

h-50.2882  MHz 

h-25.0020  MHz 

h- 84.7 185  MHz 

B  P3/2 

-h-0.4934  MHz 

h-25. 790(9)  MHz 

h- 12.4965  MHz 

CP3/2 

h-0.0560  MHz 

- 

- 

This  completes  the  section  on  the  internal  structure  of  atom  important  to  this 
paper.  Understanding  the  simple  Hydrogen  system  allowed  us  to  adapt  equations  to  fit  the 
alkali  elements.  From  a  coarse  perspective  we  saw  that  the  alkali  metal  energy  levels  can 
be  adapted  to  that  of  what  the  hydrogen’s  exact  levels  are.  Since  the  alkali  metals  and 
hydrogen  have  only  one  valence  electron  orbiting  in  a  ground  state  s-orbital  this  allowed 
correction  factors  to  be  inserted  to  describe  the  alkalis.  At  a  fine  level  the  intrinsic  spin  of 
the  electron  coupled  to  the  magnetic  field  created  by  the  orbiting  nucleus  (in  the 
electron’s  rest  frame).  Using  the  ideas  from  special  relativity  and  combining  them  with 
the  ideas  from  quantum  mechanics  both  spin  and  SO  coupling  were  revealed  to  be  a 
natural  phenomena  of  a  more  accurate  description  of  reality  rather  than  just  using  ad  hoc 
semi-classical  explanations.  Lastly,  the  hyperfine  structure  was  the  result  of  including  the 
magnetic  potential  into  the  Hamiltonian  of  the  electron.  In  the  case  of  hydrogen  it  was 
enough  to  use  only  the  dipole  moment  since  it  is  the  only  moment  a  proton  has.  However 
when  moving  to  the  alkali  atoms  the  nuclei  no  longer  present  a  simple  electric  and 


magnetic  field.  To  rectify  the  situation  one  has  to  begin  by  expanding  both  the  vector  and 
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electric  potentials  in  to  their  multi-pole  form.  Then  following  the  same  analysis  as  in  the 
hyperfine  hydrogen  analysis,  a  hyperfine  structure  is  revealed  to  involve  the  magnetic 
dipole  and  octupole  as  well  as  the  electric  quadrupole.  Figures  1-6  capture  the  internal 
structure  of  Cs,  Rb85,  and  Rb87  and  their  energy  levels. 
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Figure  1.  Cs  D1  Energy  Level  Diagram 
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Figure  2.  Cs  D2  Energy  Level  Diagram 
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Figure  3.  Rb85  D1  Energy  Level  Diagram 
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Figure  4.  Rb85  D2  Energy  Level  Diagram 


24 


1/2 


306.246(11)  MHz 


510.410(19)  MHz 

I 


9f  =  1/6 

(0.23  MHz/G) 


F=  2 


816.656(30)  MHz 


9f  =  -  1/6 
(0.23  MHz/G) 


F=  1 


794.978  850  9(8)  nm 
377.107  463  5(4)  THz 
12  578.950  985(13)  cm' 
1.559  590  99(6)  eV 


T 


(0.70  MHz/G) 


5l/2 


2.563  005  979  089  1 1(4)  GHz 


T 

▲ 

6.834  682  610  904  29(9)  GHz 


4.271  676  631  815  19(6)  GHz 


L 


(-0.70  MHz/G) 


Figure  5.  Rb87  D1  Energy  Level  Diagram 
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Figure  6.  Rb87  D2  Energy  Level  Diagram 
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Interaction  of  Light  and  Matter 

Time  Dependent  Perturbation  Theory  (TDPT) 

Once  the  energy  levels  in  the  alkali  atom  structures  have  been  defined  to  the 
hyperfine  level  a  number  of  questions  arise.  The  most  relevant  question  is  how  does  one 
physically  promote  an  electron  to  a  higher  energy  state?  The  obvious  answer  is  to  add 
energy  to  the  system.  We  could  collisionally  excite  the  atoms  by  bombarding  them  with 
an  inert  gas  or  perhaps  heat  the  system  up  so  that  elevated  states  become  the  new 
equilibrium.  For  the  case  of  a  DPAL  system  the  alkali  atoms  are  promoted  through  what 
is  known  as  stimulated  absorption.  We  will  consider  an  electromagnetic  wave  impingent 
on  the  alkali  atom.  A  dipole  moment  can  be  specified  between  the  valence  electron  and 
the  nucleus  and  it  is  this  dipole  moment  that  will  interact  with  electromagnetic  wave.  In 
the  case  of  a  DPAL  system  this  light  is  provided  by  the  diode  laser  pump  source.  There  is 
more  to  the  story  than  just  promotion.  Two  other  questions  naturally  arise  as  well  1.)  are 
all  hyperfine  transitions  allowed  and  2.)  how  strong  do  the  transitions  couple  to  the 
electromagnetic  radiation?  All  these  questions  can  be  answered  by  moving  to  a  time- 
dependent  regime.  This  is  necessary  because  the  system  is  no  longer  in  a  stationary  state. 
Before  the  diode  is  turned  on  the  system  is  unperturbed.  As  soon  as  the  diode  is  turned  on 
at  some  time,  t,  a  perturbing  Hamiltonian  must  be  introduced  into  the  system  for  a  length 
of  time  equal  to  how  long  the  perturbations  lasts.  Once  the  perturbation  is  gone  the 
system  will  relax  back  into  its  unperturbed  state.  TDPT  will  be  used  to  explore  the 
answers  to  these  questions. 

The  starting  point  for  TDPT  is  the  time  dependent  Schrodinger  equation.  Assume 
that  the  time  dependency  is  parameterized  by  a  strength  factor  A : 


27 


(1.35) 

h(0  =  h0+Ah1(0 

(1.36) 

The  most  important  question  that  can  be  asked  in  this  regime  is  what  is  the  probability  of 
the  initial  state,  Tj ,  transitioning  to  a  final  state,  TV?  In  the  case  of  the  excitation  of  an 
alkali  atom  what  is  the  probability  that  a  photon  will  excite  an  electron  from  the  ground 
state  to  an  excited  state? 

Start  by  expressing  the  state  ket  as  an  expansion  in  its  own  complete  basis: 

ywHiy.wK)  (1.37) 

n 

If  Hi  is  not  present  we  can  show  with  the  help  of  the  propagator  the  solution  to  cn(t)  as: 

J£cn(t)\n°)  =  U(t)\ys(0)) 

X ' Cn  (0  |  n°  )  =  X  e~'E",n'  |"°)("°  |C»(°)|"°) 
n  n 

s  j\  (  (W  —iE®t/h 

c„(0  =  c„(0)e 

However  when  Hi  is  present  there  must  be  another  term  in  the  state  function  that  evolves 
in  time.  Taking  this  into  account  the  state  function  can  be  written  as: 

\V(t))  =  'ZK(f)e-<"‘  |«")  (1.39) 
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If  the  bn(t)  change  it  is  because  of  Hi  and  not  Ho.  Equation  1.39  can  now  be  substituted 
into  equation  1.35  to  get  the  following: 


2>6„(0-.JH,(04„(0 


-iE°t/h 

e  " 


(1.40) 


Next  take  the  dot  product  of  equation  1.40  with  the  final  state  \f(t))  =  e  'E“"U  |/°^  and 
solve  for  bp 


-iE°„tlh 

e  " 


(1.41) 


Mf  =  X(/°  ^(Oln0)^^)^  (1-42) 

n 


E°  —E° 

where  (O,  =— - — . 

fn  fi 

No  approximations  up  till  now  have  been  used.  As  with  many  methods  in  quantum 
mechanics  we  want  to  see  if  higher  order  corrections  affect  zeroth  order  solutions  and  so 
the  expansion  coefficients  themselves  will  be  expanded  in  a  power  series.  This  is  the  first 
introduced  approximation  for  TDPT: 


4(0=4°’(0+4l,('M+ftf’('W+... 


(1.43) 


Once  equation  1.43  is  substituted  into  equation  1.42  we  find  that: 


|>^>«=X(/0kM 


r= 0 


0\ 

n)e  1 


oo 


r=0 


(1.44) 


Now  like  powers  of  r  can  be  equated  together  and  two  conditions  can  be  examined  with 
the  second  leading  to  a  recursion  relation. 


when  r  =  0, 


(1.45) 


when  r  ^  0, 


/m«(o=X(/°  “i(0 


(1.46) 


To  compute  the  first  order  result  the  zeroth  order  result  must  be  computed  first.  Equation 
1.45  says  that  the  zeroth  order  solution  is  a  constant  in  time  so  pick  to  solution  at  t=0.  In 
other  words  the  initial  state  before  the  perturbation  occurs: 


(1.47) 


Moving  on  to  the  first  order  result: 
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(1.48) 


^/)(0  =  X(/°  H,(0  nQ)etafibf{ 0) 


K/’HjCO 

n 

/“|h,(0 


i°\e,aft 


Integrating  equation  1.48  the  solution  for  the  time  evolved  state  up  to  first  order  is: 


bf(t)  =  b°M)+blJt) 


'fVJ  1  Vf 

-  I 

l 

h 


=  5/-ij{/»|H,(0 


o  \  "/I 


0 


dt' 


(1.49) 


The  probability  for  transition  from  initial  to  final  state  is  simply  the  absolute  value  of 
equation  1.49  squared.  Now  that  a  general  formula  for  transition  probability  has  been 
derived  the  form  of  Hi  will  make  it  applicable  to  systems  we  are  interested  in,  mainly  the 
hydrogen  atom  and  the  excitation  of  alkali  atoms  in  the  presence  of  an  electromagnetic 
field. 

The  Hamiltonian:  Zeroth  and  First  Order 

We  will  look  at  the  case  of  monochromatic  light  interacting  with  the  simple 
hydrogen  system.  Although  the  hydrogen  atom  is  neutral  it  contains  a  dipole  moment  by 
virtue  of  proton  electron  separation.  As  with  any  dipole  phenomenon  it  will  have  a 
natural  resonating  frequency  and  it  is  this  frequency  that  can  couple  to  monochromatic 
light  if  and  only  if  the  natural  frequency  and  frequency  of  radiation  match.  To  begin  we 
start  with  the  relativistically  invariant  form  of  hydrogen  in  an  electromagnetic  field: 
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+  (F(r  )  +  &/>) 


(1.50) 
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--A(r) 

c 


Y 
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where  A  is  the  vector  potential  and  0  the  electric  potential  of  the  radiation  field. 
Expanding  equation  1.50  and  substituting  into  the  Schrodinger  Equation  we  find: 


ifi  —  \y/)  -  • 
dt' ' 


2  m 


fiV  +ihV-  A  +  2ih-A-V  +  —  A' 


+  V(f)  +  0p\\yr)  (1.51) 


Utilizing  the  Coulomb  Gauge  equation  1.51  can  be  put  in  the  form  of: 


H|r}=[w„+«,(0]|r> 
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+  V(f) 
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(1.52) 


So  long  as  the  magnetic  field  is  not  absurdly  high  the  fourth  term  in  equation  1.52  can  be 
neglected.  The  form  of  the  vector  potential  is  both  oscillatory  in  space  and  in  time  and 
can  be  written  as  a  plane  wave: 

A(r,0  =  A/-“-kf)  (1.53) 


and  the  electric  field  is, 


E(r,0 


-IdA 

c  dt 


»A0e(-“-kf) 


(1.54) 
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Periodic  Perturbations 


First,  we  shall  treat  the  time  dependence  of  the  perturbation  in  equation  1.52 
separately  expressing  Hi  as: 


Rl(t)  =  Rle~cot 


(1.55) 


This  will  make  the  math  more  simplistic  since  sines  and  cosines  which  typically  describe 
wave  behavior  in  time  can  be  expressed  as  exponentials.  Suppose  at  t  =  0  the  hydrogen 
atom  comes  into  contact  with  the  perturbation  described  by  1.54  for  a  finite  amount  of 
time,  t.  This  can  now  be  handled  by  equation  1.49: 


The  probability  of  transition  is  the  absolute  value  of  equation  1.56  squared: 


(1.56) 


sin  {o)fi-toyi2 


(o)f-o))t!2 


(1.57) 


Equation  1.57  has  a  functional  form  of  sin2(x)/  x2  and  if  this  is  plotted  we  see  that  the 
probability  will  favor  a  transition  based  on  the  width  of  this  function. 
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3  2 


2  3 


The  function  has  a  full  width  half  maximum  of  about  n.  Let’s  consider  this  FWHM  value 
to  be  the  favorable  transition  threshold  so  that  we  can  relate  the  denominator  of  1.57  as: 


(1.58) 


(COf-Q))t/  2  <  K 


Multiply  equation  1.58  by  h  so  that  energy  units  are  used  and  solve  for  (On  to  yield: 


(1.59) 


As  long  as  cot  »  2n  and  we  know  how  to  calculate  initial  and  final  energy  levels  we 
only  need  to  use  a  frequency  proportional  to  the  energy  difference  to  excite  an  electron. 

Dipole  Transition 

Now  that  the  temporal  portion  of  the  perturbation  has  been  analyzed  next  the 
spatial  portions  which  are  contained  in  the  matrix  elements  of  equation  1.56  will  yield 
information  about  what  hyperfine  transitions  are  allowed  and  their  strength.  The  matrix 
elements  are  best  resolved  one  component  at  a  time. 
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(1.60) 
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Next  we  introduce  the  dipole  approximation  which  is  first  to  expand  the  exponential  in 
equation  1.60: 


Jkz 


=  \  +  ik(z-z0)-—(z-z0)2+ ... 


(1.61) 


If  the  product  of  kz  is  much  much  less  than  1  we  only  need  to  take  the  first  term  for  the 

2  ji 

exponential.  The  constant  k  is  equal  to  —  so  for  any  wavelength  that  is  greater  than 

X 

ultraviolet  this  approximation  will  remain  valid.  A  more  physical  picture  is  this:  over  the 
length  of  the  molecule  or  atom,  which  is  on  the  order  of  angstroms,  the  spatial  deviation 
of  the  electromagnetic  wave  hardly  changes  and  can  be  approximated  as  a  constant  field. 
With  this  approximation  in  mind  equation  1.60  becomes  (see  Appendix  A): 


ieh  ikz  d 
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(1.62) 


Plugging  this  into  equation  1.56  the  final  form  for  the  probability  transition  is  nearly 
complete: 
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(1.63) 


The  most  important  term  in  this  transition  probability  lies  in  the  matrix  element  of  the 
transition.  From  these  elements  we  can  understand  not  only  what  transitions  are  allowed 
between  arbitrary  initial  and  final  states,  but  also  if  the  transition  is  allowed  how  strong 
will  it  be.  The  generalized  theorem  to  calculate  these  matrix  elements  is  the  Wigner- 
Eckart  Theorem  which  states: 


™  •  I  rpq  I  .  \  {p^l  I  1  ■  \ 

a 2’J2’m2\Tk\ai’Jvm\}  = 1 - ft-T-j-y - L{j2,m2\k,q-,j„m ,) 


(1.64) 


where  the  are  spherical  tensor  operators,  the  a’s  represent  all  non-angular  quantum 

numbers,  and  j,m  are  standard  quantum  numbers.  We  see  that  the  angular  dependence  of 
the  matrix  elements  are  separated  into  the  familiar  Clebsch-Gordon  Coefficients.  This 
will  be  related  to  the  strength  of  the  transition  and  can  easily  be  calculated  from  CG 
Tables.  The  remaining  matrix  element  will  give  us  the  selection  rule.  Essentially  it  must 
be  that  in  the  reduced  matrix  element  on  right  hand  side  of  equation  1.64  j2  =  /,  ±  k  or 

else  no  transition  will  take  place.  If  we  recognize  the  following: 

Tr'-T, 1 
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(1.65) 


36 


and  substitute  r  =  x+y+z  the  Wigner-Eckart  theorem  can  be  used  to  compute  transition 


strengths  and  allowed  transitions. 

Three  Level  System 


Figure  7.  General  Three  Level  System  Diagram 


Suppose  we  have  a  three  level  atomic  system,  as  in  Figure  7,  where  the  first  and 
third  and  the  first  and  second  levels  are  coupled  by  a  dipole  allowed  transition.  Next,  the 
question  can  be  asked  what  are  the  transition  rates  associated  with  the  allowed  transitions 
from  TDPT.  Incorporated  in  these  transition  rates  are  optical  constants  known  as  the 
Einstein  A,  B  coefficients  that  have  both  classical  and  quantum  mechanical 
interpretations. 
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Einstein  A,B  Coefficients 

In  order  for  a  transition  to  take  place  from  the  first  to  the  third  level  (or  the  first  to 
the  second)  a  radiation  field  must  be  present.  Obviously  this  radiation  field  will  not  be 
composed  of  one  solitary  frequency  matched  to  the  energy  difference  of  transition,  but 
rather  will  have  a  distribution  of  frequencies  which  can  be  described  by  a  Plank  radiation 
distribution  : 


P(0J)  = 


CO 

„2  3 

K  C 


Pico 


ghco/kT  | 


(1.66) 


The  rate  at  which  atoms  in  the  first  level  are  promoted  to  the  third  level  can  be  written 
down  as: 


■^13  NM®) 


(1.67) 


where  Ni  is  the  number  density  of  atoms  in  the  first  energy  level  and  B 13  is  known  as  the 
Einstein  B  coefficient  which  has  units  of  miicol  Et .  Absorption  isn’t  the  only  process 
occurring  however.  An  excited  atom  in  the  third  state  can  do  one  of  two  things. 
Spontaneously  emit  its  absorbed  energy  in  any  direction  or  it  can  emit  stimulated 
radiation  by  means  of  a  dipole  relaxation  transition.  This  allows  us  to  write  down  a  rate 
equation  for  transition  from  the  third  level  to  the  first: 


■^31  2V3  B3lp(  co )  +  N-,  A-, ! 


(1.68) 
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where  the  A  coefficient  is  in  units  of  inverse  seconds  and  a  property  of  the  atom  itself. 
When  the  system  comes  to  equilibrium  we  can  take  advantage  of  the  Boltzmann 
distribution  of  population  levels  to  write  down  the  following  relation: 


-^3  _  Sj_ .  e-ha/kT  _  Bup((P) 

Nx  g\  -^3lP(d^)  +  ^3] 


Solving  for  the  radiation  field  we  find  a  relation  between  the  Einstein  coefficients  and  the 


electromagnetic  field  present: 
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(1.70) 


So  the  only  way  this  can  be  true  is  if  the  A,B  coefficients  have  the  following 
relationships: 
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(1.71) 


The  B  coefficient  can  also  be  understood  from  a  quantum  mechanical  point  of 
view  as  well  by  starting  with  equation  1.56.  Once  again  we  wait  a  “long  time”  so  that  the 
system  can  come  to  equilibrium  and  write  the  TDPT  coefficient  as: 
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Squaring  this  quantity  we  get  the  transition  probability: 


P^f  = 


4/r 
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S((of-co)8(cof-co) 


(1.73) 


The  double  delta  function  encountered  can  be  handled  in  the  following  way: 


88  =  lim  8(cof-  co )— !—  [  e'(a>fi  a)' dt 
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(1.74) 


Now  all  that  is  left  to  do  to  derive  the  average  transition  rate  is  to  substitute  1.74  into  1.73 
and  divide  by  T: 
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S(E°f  -E°  -hv) 


(1.75) 


If  our  initial  and  final  states  correspond  to  the  first  and  third  level  of  our  three  level 
system  and  recognizing  that  in  equation  1.67  the  quantity  Bup(co)  is  also  a  rate  we  see 
that  the  B  coefficient  expressed  through  TDPT  is: 
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The  same  analysis  can  be  done  for  the  dipole  transition  between  states  1  and  2  which  will 
be  fully  expounded  later  in  the  paper  for  the  full  rate  equation  analysis. 

Kinetic  Considerations 

Stimulated  emission  and  absorption  along  with  spontaneous  emission  are  not  the 
only  process  occurring  in  a  three  level  system.  If  there  exist  foreign  species  in  the  system 
in  a  gaseous  state,  for  instance  hydrogen  or  methane,  they  will  impart  energy  upon  the 
atomic  species  in  the  form  of  collisions.  What  these  collisions  do  to  the  states  of  the  atom 
depend  chiefly  on  two  things,  the  energy  of  the  colliding  species  and  how  well  the 
foreign  species  can  absorb  energy  in  terms  of  rotational  and  vibrational  excitations.  How 
well  a  species  can  couple  to  the  atomic  species  is  beyond  the  scope  of  this  research 
however  the  results  are  of  utmost  importance.  Suffice  to  say  a  foreign  species  can  do  the 
following: 

1)  Promote  an  atom  from  state  1  to  state  2  or  3. 

2)  Demote  and  excited  atom  from  3  or  2  down  to  1. 

3)  Regulate  the  excited  population  states. 

Each  species  will  have  its  own  ‘k’  coefficient  and  its  affect  on  the  population  levels  will 
depend  on  how  much  foreign  species  is  present  in  the  form  of  partial  pressure,  M.  These 
kinetic  factors  will  affect  each  population  rate  by  the  following  functional  form: 
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where,  x  represents  a  foreign  species  and  i,j  represents  different  energy  levels.  Figure  1 


depicts  the  three  sets  of  kinetic  interactions.  The  most  important  set  is  what  are  deemed 
the  spin-orbit  coefficients,  k23  and  k32.  These  coefficients  regulate  how  quickly 
population  is  de-excited  from  the  third  to  second  level.  Table  4  lists  the  spin  orbit 
coefficient  k32  of  Cesium  and  Rubidium  between  both  Helium  and  Ethane. 


Table  4.  Rubidium  and  Cesium  Spin  Orbit  k32  Constants 


Collision  Partners 

Value  (s'1-torr‘1) 

Cs-He 

157.601 

Cs-Ethane 

9.061-106 

Rb-He 

40137.500 

Rb-Ethane 

1.082-107 

Absorption  Cross  Section 

In  equation  1.75  the  delta  function  dictates  that  a  transition  will  only  take  place  at 
a  frequency  equal  to  the  difference  of  the  transition  energy  levels.  On  the  outset  it  seems 
that  only  one  frequency  will  promote  a  transition,  but  this  is  not  entirely  true.  Just  as  in 
blackbody  radiation  we  find  a  spectrum  of  frequencies  peaked  at  a  certain  value,  there  is 
also  a  spectrum  of  frequencies  that  can  be  absorbed  by  a  stimulated  response  peaked  at 
E f  -  Ej .  What  modifies  this  absorption  profile  are  deviations  of  an  atom  from  its 

quantum  mechanically  defined  energy  level.  These  deviations  come  in  the  form  of 
natural,  thermal,  and  pressure  deviations  that  tend  to  broaden  the  frequencies  of  radiation 
which  would  promote  an  atom. 

The  defining  equation  that  entails  all  of  this  information  will  be  termed  the 
differential  absorption  law: 
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(1.78) 


1  dl(CD) 
I  (CD)  dx 


-  -AN  ■  c0  •  g(co) 


where  AN  is  the  population  inversion  of  the  absorption  manifold.  Typically  given  in 
units  of  cm"  is  the  integrated  cross  section  c0  and  is  modified  in  frequency  space  by 
g(co )  or  the  lineshape  of  absorption.  The  lineshape  function  will  capture  all  the  effects 
that  broaden  the  absorption  profile. 

Integrated  Cross  Section 

The  cross  section  can  be  thought  of  as  a  measure  proportional  to  how  much 
radiation  is  absorbed  by  an  ensemble  of  atoms  at  a  particular  frequency.  A  relation  can  be 
made  between  the  B  coefficient  and  the  cross  section.  The  total  rate  at  which  photons  are 
absorbed  due  to  a  stimulated  response  in  the  range  of  co  to  to  +  dro  is: 


Rudco  =  Bng{co)p((o){Nx  -  —  N3)do)  (1-79) 

S  3 

Multiplying  this  rate  by  ti(0  and  supposing  that  the  power  density  is  due  to  a  beam  of 
cross-sectional  area  A  we  find  the  power  lost  due  to  absorption  as: 


-A P  =  hcoBug(co)p(co)(Nl-—N3)AdcoAx  (1.80) 

ft 


where  Ax  is  simply  how  far  the  beam  propagates.  Dividing  both  sides  of  equation  1.80 
and  noting  that  I(co)  =  p(co)c  we  find  that: 
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(1.81) 


Comparing  equations  1.81  with  the  differential  law  of  absorption,  equation  1.78,  and 
using  the  relationship  between  the  A,  B  coefficients  in  equation  1.71  we  can  calculate  the 
integrated  cross  section  from  simply  knowing  the  radiative  lifetime: 


=  ^M 
®  & 


(1.82) 


Natural  Broadening 

Let  us  assume  that  for  a  moment  we  only  have  a  two  level  system  with  one  atom 
in  the  upper  level  at  a  strictly  defined  energy.  According  to  the  uncertainty  principle  if 
this  atom  is  indeed  so  energetically  well  defined  the  time  it  would  exist  in  that  state 
would  be  nearly  infinite!  But  the  fact  of  the  matter  is  that  atoms  do  radiate  and  so  it  must 
be  that  our  energy  levels  aren’t  defined  to  the  precision  we  have  made  them  out  to  be. 
Luckily  this  has  not  affected  our  analysis  to  this  point,  but  to  proceed  further  we  can  no 
longer  ignore  Heisenberg’s  Uncertainty  Principle.  Tightly  clustered  around  the  upper 
level  is  a  small  spectrum  of  energy  levels  that  an  excited  atom  can  exist  in. 
Correspondingly  there  is  a  spectrum  of  energy  levels  about  the  lower  state  as  well. 
Pictorially  we  can  have  a  transition  from  the  lower  band  of  energies  to  the  upper  band  of 
energies: 
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Figure  8.  Effect  of  Uncertainty  Principle  on  Energy  Levels 


Mathematically  this  distribution  of  absorption  frequencies  is  Lorentzian  in  nature: 


1 


gO-V2l)=^ 


Sv 


2k  (v-v2l)2  +(8v/  2)2 


(1.83) 


where  c)vis  the  Full  Width  Half  Maximum  (FWHM)  and  is  equal  to  the  A  coefficient 
divided  by  2%. 

Doppler  Broadening 

Another  phenomenon  that  will  cause  the  absorption  spectra  to  broaden  as  well  is 
the  Doppler  effect.  At  a  certain  temperature  the  atoms  to  be  excited  will  move  with  a 
certain  kinetic  energy  dictated  by  a  Maxwell  Boltzmann  distribution.  Some  fraction  of 
the  atoms  in  the  system  will  therefore  be  moving  parallel,  anti-parallel,  or  some 
combination  to  the  incoming  electromagnetic  wave.  As  such  in  the  frame  of  the  moving 
atom  the  imparted  EM  wave  will  be  red  or  blue  shifted  accordingly.  This  shifting  must 
then  match  an  absorbing  eigenfrequency  if  the  atom  is  to  be  promoted.  This  leads  to  a 
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Gaussian  distribution  in  the  absorbing  lineshape  of  the  atoms  with  a  corresponding 
FWHM  dependent  on  temperature: 


g  (v-v21) 


_2_  / 1  n  ( 2  ) 
Sv\  n 


,e-41n(2  )[(v-v21)/5v]2 


t)v  —  2v2I 


2fc?Tn(2) 

me1 


(1.84) 


Pressure  Broadening  and  Pressure  Shifting 

The  last  observable  fact  that  we  will  consider  is  pressure  broadening,  sometimes 
termed  collisional  broadening,  and  pressure  shifting.  A  pure  quantum  mechanical 
treatment  is  very  difficult  to  obtain  for  it  depends  on  interatomic/intermolecular 
potentials  between  colliding  parties.  When  two  collisional  partners  are  far  apart  their 
energy  levels  are  that  of  their  unperturbed  state.  However  as  they  come  within  collision 
distances  the  potential  surfaces  of  both  the  lower  and  upper  level  begin  to  change.  The 
absorption  frequency  then  becomes  a  function  of  intermolecular  distance.  The  higher  the 
pressure  the  more  the  atomic  potential  surfaces  are  effected  which  serves  to  broaden  the 
frequency  of  absorption.  Correspondingly  high  pressures  can  also  totally  shift  where  the 
transition  frequency  occurs  by  shifting  both  the  upper  and  lower  level  unperturbed  energy 
levels.  Both  the  broadening  coefficient  and  shifting  coefficient  depend  on  the  particular 
collision  partner.  No  closed  form  quantum  mechanical  derivation  exists  for  such 
complicated  potentials  as  such  these  parameters  are  experimentally  measured  and 
reported.  For  many  cases  it  has  been  experimentally  verified  that  a  Lorentzian  profile  best 
approximates  this  effect: 


46 


(1.85) 


21  2k  (v  — (v21  +  Av))2  +  (<5v/2)2 

<5v  =  X^x 

A  y='Lb',P, 

X 

where  b  is  the  broadening  coefficient,  b’  is  the  shifting  coefficient,  and  P  is  the  pressure 
in  torr  of  each  separate  collision  species. 

Complete  Set  of  Rate  Equations 

Now  one  can  write  down  the  full  dynamics  of  a  three  level  atomic  system  in  the 
presence  of  a  stimulating  radiation  field  and  foreign  collision  partners  as: 


1*  =  Ij^hTL  *  N,B"y ~  NAl  -  N' '  K'-M- + £  *”■ M’ ) 


+N2  ■  ^  &23,  M x  ^  k[3^  Mx 

X  X 

^L  =  L^2L(-Nibx+NA2)-N1A1i -N2  -(X^,  ,M,  +  £**  ,M,)  (1.86) 


+a v5X-*/.+jv2Xm- 


J7V, 


i  _ 


f  c/M  diV,  ^ 
— -  +  — - 
dt  dt 


Saturation  Intensity 

The  saturation  intensity  is  defined  to  be  the  intensity  that  decreases  the  absorption 
or  gain  to  exactly  half  that  of  the  small  signal  gain.  The  value  of  this  intensity  is 
particular  for  each  system  at  specific  system  parameters  for  it  depends  on  the  cross 
section  which  is  a  function  of  frequency,  temperature,  and  pressure. 
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(1.87) 


O(C0) 

In  fact  equation  1.78  is  a  special  case,  known  as  the  Beer’s  Law  Limit,  of  the  general 
differential  absorption  law.  The  general  differential  absorption  law  contains  the  saturation 
intensity  as  follows: 


1  dHt»)_-N-a0-g(a>)-a  , 

Ha>)  dx  p  +  x!M  '  ' 

s 

with  functions  a,  ft,  and  X  determined  from  the  specific  system  in  question  (See 
Appendix  A).  When  the  input  intensity  is  far  below  the  saturation  intensity  we  recover 
equation  1.78.  If  the  input  intensity  is  far  above  the  saturation  intensity  then  the  first  term 
in  the  denominator  of  1.88  can  be  ignored  and  we  find  linear  behavior: 


dl(co)  _  -N  ■  Is  ■  <Jn  •  g(co)  ■  a 
dx  x 


(1.89) 


Summary 

With  this  background  a  comprehensive  model  can  be  pieced  together  to  explore 
the  various  facets  of  what  would  affect  a  three  level  system  when  it  is  absorbing  and 
lasing.  The  task  at  hand  is  to  integrate  all  these  separate  phenomena  and  then  check 
against  experimental  results  to  see  the  accuracy  of  the  model.  Once  this  is  performed  how 
important  is  one  parameter  over  the  other,  or  in  other  words  how  sensitive  is  the  system 
to  certain  parameters,  can  be  explored  in  the  model  rather  than  in  experimental 
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laboratories.  This  model  will  also  serve  as  a  tool  to  help  fine  tune  what  a  final  DPAL 
system  might  look  like. 
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II.  Description  of  DPAL  Absorption  Model 


Chapter  Overview 

This  chapter  describes  how  the  absorption  and  lasing  model  of  a  three  level 
DPAL  system  is  constructed  and  controlled.  The  two  alkalis  modeled  in  this  system  are 
Rubidium  and  Cesium  with  the  energy  level  diagram  as  seen  below: 


This  model  is  completely  written  using  Wolfram’s  Mathematica  6.0. 1.0  32-bit 
version.  This  is  a  totally  symbolic  math  package  capable  of  solving  and  plotting  user 
defined  equations.  The  model  first  begins  by  the  user  defining  the  parameters  of  the 
system:  temperature  of  the  absorbing  cell,  temperature  of  the  alkali  cell,  pressure  of  the 
foreign  species  (helium  and  ethane),  power  of  the  pump  source,  and  finally  the  full  width 
half-maximum  of  the  pump  source  modeled  as  either  a  Gaussian  or  square  wave.  One 
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simply  uses  the  control  panel  for  absorption  which  consists  of  a  set  of  “sliders”  to  pick 
what  value  one  wishes  for  each  parameter.  Below  is  a  picture  of  what  the  temperature 
parameter  input  sliders  looks  like. 

■  Control  Panel  of  parameters  for  Absorption  : 

■  Enter  Temperature  of  Cell  (K)  : 

{Slider  [Dynamic [Tl]  ,  (0,  1000,  1},  ImageSizc  -*  500,  ContinuousAction  -♦  False] ,  Dynamic [T1 ] } 

{■_  ..  =D-  -  242} 

■  Temperature  of  Cesium  Sample  (K) : 

{Slider  [Dynamic  [T2]  ,  {200,  500,  1},  Image  Size  -»  500]  ,  Dynamic  [T2  ]  ) 

{O  200} 

*  note  the  melting  point  for  Cs  is  301 .59  K  * 

Figure  10.  Parameter  Control  Panel 

Once  the  parameters  are  set  all  that  needs  to  be  done  is  run  the  Mathematica  kernel  and 
the  results  are  output  at  the  end  of  the  program  code. 

Constants  and  properties  of  the  System 

All  physical  constants  are  taken  from  the  Committee  on  Data  for  Science  and 
Technology’s  (CODATA)  recommendation  list  for  fundamental  physical  constants  dated 
2006.  The  physical  properties  and  optical  constants  of  Rubidium  and  Cesium  are  taken 
from  the  experimental  work  of  David  A.  Steck  at  the  Los  Alamos  National  Laboratory. 

Computing  Energy  Levels  and  Transitions: 

The  first  step  to  computing  what  frequencies  will  induce  a  transition  is  to  first 
compute  the  available  energy  levels  from  quantum  mechanics.  Equations  1.4,  1.22,  and 
1.34  allow  the  calculation  of  energy  levels  to  the  hyperfine  structure.  Designate  the 
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transition  from  the  ground  state,  6  2S]n ,  to  the  excited  state  ,  6 2 Pu2 ,  the  D1  transition 

manifold  and  from  the  ground  state  to  the  next  higher  excited  state,  6  2P3/2 ,  the  D2 

transition  manifold.  Using  the  rules  of  the  addition  of  angular  momentum,  in  the  case  of 
Cesium,  the  ground  state  F  value,  where  F  equals  the  addition  of  electron  spin,  angular 
momentum,  and  nuclear  spin,  takes  on  3  or  4.  Then  understanding  equation  1.64,  that  is 
AF  =  0,  ±1,  we  can  calculate  at  what  frequency  to  expect  a  hyperfine  transition  in  both 
the  D1  and  D2  manifolds.  Tables  5-7  predict  these  hyperfine  transition  frequencies.  Note 
that  in  Tables  5-7  the  transition  frequency  of  the  fine  structure  transitions  or  line  center, 

S 1/2  — »  P 1/2  and  S1/2  — ►  P3/2,  have  been  offset  to  zero.  This  is  why  there  appears  to  be  a 
negative  transition  frequencies. 


Table  5.  Cesium  Transition  Frequencies 


D1  Manifold 

D2  Manifold 

Transition 
(Gnd  — >  Excited) 

Frequency 

(MHz) 

Transition 
(Gnd  — » Excited) 

Frequency 

(MHz) 

3— »3 

4514.04 

3— >2 

4829.68 

3— >4 

5681.72 

3-4 

4984.81 

4-4 

-  4678.60 

3-4 

5182.19 

4— >4 

-3510.92 

4-4 

-  4207.82 

4—4- 

-  4010.44 

4— >5 

-  3757.58 

Table  6.  Rubidium  85  Transition  Frequencies 


D1  Manifold 

D2  Manifold 

Transition 
(Gnd  — *  Excited) 

Frequency 

(MHz) 

Transition 
(Gnd  — >  Excited) 

Frequency 

(MHz) 

2-4 

1559.92 

2—4 

1657.64 

2-4 

1921.5 

2— >2 

1687.01 

3-4 

-1475.81 

2— >3 

1750.41 

3— >3 

-1114.23 

3-4 

-1348.72 

3-4 

-1285.32 

3-4 

-1164.68 
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Table  7.  Rubidium  87  Transition  Frequencies 


D1  Manifold 

D2  Manifold 

Transition 
(Gnd  — *  Excited) 

Frequency 

(MHz) 

Transition 
(Gnd  — *  Excited) 

Frequency 

(MHz) 

1->1 

3761.27 

1-^0 

3969.6 

1  —>2 

4577.92 

1->1 

4041.82 

2— >1 

-3073.42 

l->2 

4198.77 

2->2 

-2256.76 

2— H 

-2792.86 

2— *2 

-2635.92 

2^3 

-2369.27 

Computing  Lineshapes 

If  the  frequencies  computed  for  the  hyperfine  transitions  of  the  D1  and  D2 
manifolds  are  plotted  we  would  produce  what  is  known  as  a  stick  spectrum.  As  discussed 
earlier  in  this  paper  we  have  to  account  for  three  phenomena  that  would  broaden  this 
stick  spectrum:  natural,  Doppler,  and  pressure.  Not  only  that  we  also  know  that  each 
transition  has  a  certain  strength  calculated  from  the  matrix  elements  via  TDPT.  In  order 
to  do  this  we  form  what  is  known  as  the  Voigt  profile.  This  is  simply  the  convolution  of  a 
Gaussian  with  Lorentzian.  The  formed  used  of  the  Voigt  profile  in  this  model  is  as 
follows: 


V(v,v,)  = 


Re[g~z(v,v,)2  -Erf  ifr.v,))] 


5vD 

V81o§(2) 


(1.90) 
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T 
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V/2\ 


,8vL 
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2 


(1.91) 
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where  the  full  width  half  max  of  the  Gaussian  and  Lorentzian  profiles  are: 


8vD  =  (vt  +vDlorD2)' 


/  \1/2 

'  SkTLog(2)  A 


SvL  = 


- h  yH  e  ■  PHe 

2k 


y  ™Cs'c 

(  J 


v  294. 15  j 


J 

A1/2  /  T 

+  y  Me  PM  A 


yl,2\ 


(1.92) 


v  294.15  j 


Each  hyperfine  transition,  Vt ,  will  be  broadened  according  to  equations  1.90-1.92.  The 
summation  of  all  transitions  in  a  particular  manifold  will  be  the  lineshape  function  g(v) . 
Each  transition  is  further  modified  by  the  strength  of  transition,  Table  8-10,  and  the 
relative  population  of  the  hyperfine  split  ground  state  based  upon  a  Boltzmann 
distribution. 


g(^  =  l§S^V(v,v,) 

t  ^ 


(1.93) 


Table  8.  Relative  Strength  of  Transition  for  Cesium 


D1  Manifold 

D2  Manifold 

Transition 
(Gnd  — *  Excited) 

Relative  Strength 

w 

Transition 
(Gnd  — >■  Excited) 

Relative  Strength 

w 

3->3 

1/4 

3— >2 

5/14 

3— >4 

3/4 

3— G 

3/8 

4->3 

7/12 

3— >4 

15/56 

4— >4 

5/12 

4->3 

1H2 

4->4 

7/24 

4— >5 

11/18 
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Table  9.  Relative  Strength  of  Transition  for  Rubidium  85 


D1  Manifold 

D2  Manifold 

Transition 
(Gnd  — *  Excited) 

Relative  Strength 

(St) 

Transition 
(Gnd  — » Excited) 

Relative  Strength 

(St) 

1-G 

1/6 

1-d) 

1/6 

1— >2 

5/6 

1-G 

5/12 

2-G 

1/2 

1  — >2 

5/12 

2— »2 

1/2 

2— G 

1/20 

2->2 

1/4 

2— *3 

7/10 

Table  10.  Relative  Strength  of  Transition  for  Rubidium  87 


D1  Manifold 

D2  Manifold 

Transition 
(Gnd  — *  Excited) 

Relative  Strength 

(St) 

Transition 
(Gnd  — *  Excited) 

Relative  Strength 

(St) 

1-G 

1/6 

1-G) 

1/6 

1— >2 

5/6 

1-G 

5/12 

2— G 

1/2 

1  — >2 

5/12 

2->2 

1/2 

2— >1 

1/20 

2-^2 

1/4 

2->3 

7/10 

Table  11.  Broadening  and  Shifting  Coefficients  for  Cesium 


D1  Manifold 

D2  Manifold 

Species 

Broadening 

Coefficient 

(MHz/torr) 

Species 

Shifting  Coefficient 
(MHz/torr) 

yHe 

26.21 

yHe 

23.50 

yMe 

29.1 

yMe 

26.2 

SHe 

4.46 

SHe 

0.75 

SMe 

-7.92 

SMe 

-7.00 

(No  experimental  values  exist  for  Cs-Ethane,  so  Rb-Methane  values  are  used  instead) 


Solving  the  Rate  Equations  and  Absorption  of  Pump  Intensity 

The  form  of  the  population  inversion,  N3-  (g3/gi)Ni,  is  needed  to  understand  how 
pump  intensity  is  absorbed  or  attenuated  through  the  absorption  cell..  Assuming  that  the 
input  laser  is  impingent  upon  the  cell  of  alkali  for  a  time  much  longer  than  any  lifetime  of 


55 


the  system  a  steady  state  regime  can  be  used  to  solve  the  rate  equations.  The  rate 
equations  then  give  how  the  population  density  numbers  of  each  level  change  with 
respect  to  system  parameters.  First,  set  equations  1.86  such  that  dN/dt=0,  i,e,  steady  state 
conditions.  Secondly,  add  the  constraint,  the  total  number  of  alkali  atoms,  N,  do  not 
change,  into  the  system  of  equations.  This  is  simply  stated  as  N2  =  N  -  N1  -  N3  and  is 

substituted  into  the  time  rate  of  change  of  population  N3.  Lastly,  rearrange  the  equations 
in  to  a  matrix/vector  form  such  that  the  “vector”  is  comprised  of  the  population  densities. 
N3 .  As  such  the  matrix  equation  becomes: 


/  1 

p  31  x 


- k. 
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hv  T, 

p  3 


\32 
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hv  T, 

p  3 


^21  1 

hvL  T2 


-k  +  7'a-  1 

K  23  + 


hv, 


I  s 
hvP 

K  ,  I l? 21  Si 


v12 


8l 


+  (k]2  +  ^13 ) 


^iPi\  8 2  83 


hvL  8i  hvP  8 1 


r-NK 3  ' 

n2 

- 

0 

W 

0 

k  7 

(1.94) 


where  the  shorthand  notation  k  -'S'  M  k.  is  used  for  the  spin  orbit  constants  (reference 

ij  x  ij  x  A 

X 

equation  1.77).  Using  Kramer’s  Rule  for  a  3x3  matrix  the  vector  in  equation  1.95  can  be 
solved  for  and  the  components  of  this  vector  are  the  quantities  sought  after  (Appendix  A). 

Absorption  Profiles 

The  absorption  profile  is  derived  from  equation  1.78  and  Mathematica  will  be 
used  to  plot  the  attenuation  of  the  input  pump  as  a  function  of  frequency  and  distance  into 
the  absorbing  cell.  The  population  inversion  used  to  solve  equation  1.78  was  derived 
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from  solving  equation  1.95.  In  computing  population  inversions  the  degeneracy  of  each 
state  must  be  factored  in  as  well.  For  the  D2  manifold  the  inversion  density  is: 


AN  -  N3-  —  Nl 
Si 


(1.95) 


Before  solving  equation  1.88  one  of  three  scenarios  can  be  chosen.  First,  absorption  and 
lasing  can  be  turned  “off’  by  setting  /  and  IL  equal  to  zero.  This  would  say  only  kinetic 


processes.  Second,  setting  IL  to  zero  would  only  allow  for  stimulated  absorption  and 

kinetic  processes  to  run  in  the  simulation.  Lastly,  by  leaving  the  rate  equations  as  is  all 
processes  are  allowed  to  run.  As  the  particular  form  of  the  solution  will  be  slightly 
different  based  on  what  processes  are  allowed  to  run,  the  attenuation  of  the  pump  beam, 
using  Mathematica’s  solver  techniques,  will  have  a  functional  form  similar  to  the 
following: 


/(x,v)  =  Isat(v)  •  ProductLog 


^  -^-—AN  x  cjj  j  (v ) 

0*sat  .  T 

e  J0 


4» 


(1.96) 


where  Io  is  the  input  power  set  by  the  user  in  the  parameter  control  sliders.  This  model 
uses  both  a  Gaussian  and  square  wave  input  profile  which  also  has  a  user  defined  full 
width  half  maximum. 

The  Product  Log,  also  known  as  the  Lambert-W  function,  seeks  to  invert  the 
equation  / ( w )  =  we"'  and  solve  for  w.  This  complex  transcendental  function  is  not 
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injective  and  therefore  will  have  more  than  one  solution.  However,  since  the  parameters 
of  the  system  are  user  defined,  it  will  be  obvious  when  a  particular  value  of  equation  1.97 
is  unphysical.  When  this  happens  Mathematica  can  be  used  to  search  for  physical 
solutions. 

Differential  Transmittance 

The  last  aspect  that  this  model  simulates  is  what  is  termed  to  be  the  differential 
transmittance.  The  differential  transmittance  is  defined  as  the  percent  chance,  that  once 
steady  state  conditions  have  been  established,  a  photon  propagated  through  the  entrance 
of  the  cell  will  be  transmitted  at  a  particular  distance  in  to  the  cell.  We  say  that  the  cell  is 
transparent  from  the  entrance  of  the  cell  to  the  distance  at  which  the  differential 
transmittance  equals  zero  percent.  Beyond  this  point  the  cell  becomes  totally  opaque.  The 
boundary  between  these  two  regions,  transparency  and  opaqueness,  is  called  the  bleached 
wave  front.  As  such  saying  a  cell  is  completely  bleached  is  equivalent  to  saying  the  cell  is 
totally  transparent.  If  a  cell  is  not  completely  bleached  then  there  exists  a  region  in  the 
cell  that  is  opaque  to  incoming  photons.  This  differential  transmittance  therefore  is  a 
visual  representation  of  what  a  bleached  wave  looks  like  and  can  be  used  to  determine 
how  far  a  bleached  wave  penetrates  a  particular  absorbing  cell.  The  functional  form  of 
differential  transmittance  is  derived  from  taking  the  derivative  of  equation  1.97  with 
respect  to  the  initial  pump  intensity. 

dl(x,v )  _  c|transmjttancc('x^v/j  (1.97) 

dl0(0,v) 
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Summary 

Using  the  principles  put  forth  in  the  beginning  of  this  paper  this  model  simulates 
how  an  initial  pump  beam  is  attenuated  through  an  absorbing  cell  of  alkali  vapor  under  a 
variety  of  adjustable  system  parameters.  Every  aspect  of  the  model,  energy  levels, 
saturation  intensities,  lifetimes,  etc.  can  be  probed  in  detail  to  help  understand  what 
parameters  the  system  is  most  sensitive  to  and  which  parameters  have  subtle  of  drastic 
effects  on  input  attenuation.  Also  with  such  a  detailed  model  analytic  solutions  derived 
by  other  means  can  be  compared  to  see  the  efficacy  of  its  assumptions. 
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III.  Model  Results 


Chapter  Overview 

This  chapter  will  focus  on  what  type  of  outputs  this  DPAL  absorption  model  can 
produce  along  with  its  comparisons  to  experimental  results  and  analytical  derivations. 
Also  we  can  show  how  various  features,  like  cross  section  and  saturation  intensities, 
change  with  change  in  system  parameters.  The  focus  of  the  results  will  center  on  two 
different  regimes.  The  first  will  be  akin  to  what  we  would  see  in  a  Beer’s  Law 
experiment:  low  power,  low  alkali  concentration,  and  a  low  amount  foreign  species  in  the 
absorption  cell.  The  second  regime  will  be  more  of  what  would  occur  in  a  typical  DPAL 
absorbing  system:  high  input  power  and  high  Helium/Ethane  pressure.  Two  general 
comparisons  will  be  made  as  well.  The  first  will  be  a  comparison  of  this  numerical  model 
with  an  analytical  quasi-two  level  analysis  (QTLA)  model.  The  second  comparison  will 
be  how  attenuation  is  affected  in  the  presence  of  a  lasing  field  solely  within  the  quasi-two 
level  regime.  Table  1 1  contains  the  parameters  that  will  be  used  for  the  two  regimes 
discussed. 


Table  12.  Simulation  Parameters 


Parameter 

Beer’s  Law  Regime 

DPAL/QTLA  Regime 

Temperature  of  Cell 
(K) 

400 

455 

Alkali  Concentration 
(#/volume) 

6.7x  10" 

6.1  x  1012 

He  Pressure 
(torr) 

50 

200 

Ethane  Pressure 
(torr) 

50 

600 

Input  Peak  Intensity 
(W/cm2) 

.60 

200 

60 


Spectroscopic  Comparison 

Step  by  step  analysis  and  comparison  of  this  model  must  start  from  the  beginning 
if  the  model  is  to  be  viable.  The  first  step  in  this  process  is  to  simply  compare 
experimental  spectroscopic  results  with  the  lineshape  of  this  model.  From  equation  1.78 
we  can  write  the  solution  as: 


Ln 


fI_' 
V7°  J 


-NL-o(y) 


(1.98) 


where  L  is  the  length  of  the  cell.  The  right  hand  side  of  equation  1.99  is  simulated  and 
compared  to  a  transmission  experiment.  The  results  are  in  Figure  1 1  below. 


Simulated  vs  Experimental  Lineshape 
I 


10000 


MH/ 


Figure  11.  Simulated  vs.  Experimental  Lineshape 


Predicted  Cross  Sections  and  Pressure  effects: 

As  discussed  in  the  background  section  there  are  three  main  processes  that  lead  to 
a  broadening  of  the  lineshape:  natural,  Doppler,  and  pressure  broadening/shifting.  Of  the 


three  pressure  broadening/shifting  has  the  most  dramatic  effect  on  line  shape  as  captured 
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in  Figures  12-17.  As  the  pressure  rises  the  hyperfine  peaks  begin  to  merge  into  one 
another.  The  cross  section  at  line  center  rises  until  it  reaches  a  maximum  at  slightly  more 
than  600  torr  of  foreign  gas  species.  At  this  point  all  hyperfine  transition  lines  have 
blended  and  with  the  added  pressure  this  value  begins  to  descend  again.  Also  that  once 
the  line  shape  has  become  one  broad  peak  that  the  pinnacle  of  this  peak  is  being  shifted 
toward  positive  frequency  in  Figures  15-17.  As  the  pressure  keeps  increasing  the  two 
trends,  broadening  and  shifting,  become  more  pronounced  and  more  dominate.  Natural 
and  Doppler  broadening  are  still  accounted  for,  however  the  pressure  effects,  even  when 
the  pressure  is  low,  tends  to  overshadow  these  two  processes.  In  practice  and  in  a  real 
DPAL  system  pressure  broadening  will  be  the  chief  mechanism  to  broaden  the  cross 
section  so  that  it  matches  the  input  shape  of  the  pump. 
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Figure  12.  031  Cross  Section  at  0  torr  He  and  0  torr  Me 
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cr 3i  Cross  Section 
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Figure  13.  031  Cross  Section  at  100  torr  He  and  100  torr  Me 
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Figure  14.  031  Cross  Section  at  200  torr  He  and  200  torr  Me 
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cr 3i  Cross  Section 
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Figure  15.  031  Cross  Section  at  300  torr  He  and  300  torr  Me 


cr31  Cross  Section 
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Figure  16.  031  Cross  Section  at  600  torr  He  and  600  torr  Me 
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cr3 1  Cross  Section 
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Figure  17.  031  Cross  Section  at  1000  torr  He  and  1000  torr  Me 


Lifetimes  and  Saturation  Intensities: 

This  model  takes  spectral  look  into  the  makeup  of  the  saturation  intensity.  As 
described  by  equation  1.87  the  spectral  saturation  intensity  is  inversely  proportional  to 
the  lineshape.  Just  as  the  lineshape  varies  dramatically  with  pressure  so  does  the  spectral 
saturation  intensity.  The  value  of  the  saturation  intensity  is  of  utmost  importance  because 
of  its  relation  to  analytical  forms  and  approximations  as  seen  in  the  general  differential 
absorption  law.  Three  methods  were  employed  to  compute  the  value  of  the  average 
saturation  intensity  from  the  spectral  saturation  intensity.  The  first  was  to  take  an  average 
over  a  user  defined  frequency  interval.  The  second  method  was  to  take  the  value  of  the 
spectral  saturation  intensity  at  each  hyperfine  transition  (six  for  the  D2  manifold)  and 
simply  average  those  values.  Lastly,  one  can  average  over  the  lineshape  (which  is  unity) 
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and  simply  use  the  value  for  the  integrated  cross  section  for  equation  1.87.  In  the  end  the 
second  method  was  used  to  compute  the  average  saturation  intensity.  The  first  method 
can  become  cumbersome  and  computationally  lengthy  if  too  big  a  frequency  interval  is 
chosen.  Also  there  is  no  clear  delineation  of  where  one  should  mark  the  beginning  and 
the  end  of  the  frequency  interval.  The  third  method  rids  the  saturation  intensity  of  any 
dependence  on  system  parameters  and  it  simply  becomes  a  constant.  Method  number  2 
seems  to  be  the  most  sensible  and  logical  for  it  has  two  advantages:  1.)  it  takes  into 
account  pressure  variations  in  the  lifetime  of  the  state  and  cross  section  and  2.)  the  six 
most  important  frequencies  are  undoubtedly  accounted  for  and  given  importance.  Figures 
18  and  19  show  plots  of  the  spectral  dependence  of  the  saturation  frequencies  in  both 
regimes  and  Table  13  displays  calculated  averages  for  all  three  methods. 
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Figure  18.  Spectral  Saturation  Intensity  for  Beer’s  Law  Regime 
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Spectral  Saturation  Intensity  of  DPAL  Regime 
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Figure.  19.  Spectral  Saturation  Intensity  for  DPAL  Regime 


Table  13.  Average  Saturation  Intensities 


Isat  (W/cm2) 

Method  #1 

Method  #2 

Method  #3 

Beer’s  Law  Regime: 

7.80 

4.05 

0.00040 

DPAL  Regime: 

24.93 

21.34 

0.00040 

Absorption  Profiles 

Rate  Equation  Population  Inversions: 

In  order  to  solve  equation  1.78  the  population  inversion  is  first  calculated.  Using 
equations  1.95  and  1.96  the  population  inversion  for  absorption  is: 


AA31=- 


N{  2  +  £|3T3  +  k 23^] 2^2 ^3  -k^k^Tfl^) 


(1.99) 


(1  +  k]2T2  +  {k[2  +  k23  (k\2  k32  )t2  +  kl3k32x2  )t3  )  +  (3  +  kl2r2  +  1k32T2  )A3  3t3 

Isat 
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Once  equation  1.101  is  substituted  into  1.78  and  all  necessary  functions  and  parameters 
are  defined  the  solution  for  absorbed  intensity  can  be  explored  in  detail  from  a  spatial  and 
spectral  perspective.  As  mentioned  before  when  the  input  intensity  is  much  greater  than 
the  saturation  intensity  linear  behavior  is  expected.  When  below  the  saturation  intensity 
exponential  behavior  is  expected.  The  full  solution  is  then  a  non-linear  mixture  of  the 
two. 

Frequency  and  spatial  representation  of  Intensity 

With  equation  1.78  solved  we  first  start  with  the  Beer’s  Law  regime  and  plot  what 
the  absorption  profile  looks  like  as  a  function  of  space  and  frequency.  A  square  wave 
input  pump  with  a  FWHM  of  20,000  MHz  was  used  with  a  peak  intensity  of  0.60 
Watts/cm".  Noting  that  the  saturation  frequency  for  this  regime  was  4.05  Watts/cm“  the 
input  is  well  below  the  saturation  value.  Figure  20  displays  the  solution  to  equation  1.78. 


3D  Representation  of  Beer's  Law  Regime  Absorption 
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Figure  20.  3D  Representation  of  Beer’s  Law  Regime  Absorption 


With  a  pretty  low  system  pressure  the  cross  section  is  highly  peaked  around  its 
hyperfine  transitions.  Also  since  the  input  is  below  the  saturation  frequency  we  would 
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expect  the  intensity  to  decay  exponentially  in  space  and  be  centered  on  the  hyperfine 
peaks.  Indeed  this  is  what  is  captured  in  Figure  20.  The  two  valleys  do  decay 
exponentially  showing  that  all  input  power  has  been  absorbed.  Upon  looking  at  Figure 
20  one  might  also  ask  why  there  is  no  decay  in  between  the  hyperfine  valleys  since  the 
saturation  frequency  is  so  high?  The  saturation  frequency  is  inversely  proportional  to  the 
cross  section.  So  a  low  cross  section  means  high  saturation  frequency.  However  in  this 
case  the  cross  section  is  so  low  between  the  valleys  that  hardly  any  absorption  takes 
place.  So  the  saturation  intensity  has  no  effect  in  this  region.  The  region  is  therefore 
completely  transparent.  An  electromagnetic  wave  with  a  bandwidth  of  around  7000  MHz 
centered  about  the  fine  transition  frequency  would  simply  transmit  through.  This  window 
has  an  interesting  feature  in  that  its  bandwidth  can  be  completely  controlled  by  the 
pressure  of  the  system.  Increasing  the  pressure  would  cause  the  bandwidth  to  decrease. 

The  plot  for  the  DPAL  regime  looks  different  than  the  Beer’s  Law  regime.  With 
the  system  pressure  being  a  total  of  800  torr  the  hyperfine  peaks  have  been  broadened  to 
the  point  where  the  cross  section  is  one  peak  akin  to  Figure  15.  The  input  pump  for  this 
simulation  is  a  Gaussian  with  a  FWHM  of  20,000  MHz  and  a  peak  intensity  of  200 
W/cnT.  Being  way  above  the  saturation  frequency  linear  behavior  is  expected.  Figures  21 
and  22  are  the  model  results. 
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Figure  21.  Space  and  Frequency  Representation  of  DPAL  Regime  Absorption  (Head  on) 


Space  and  Frequency  Representation  of  DPAL  Regime  Absorption 
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Figure  22.  Absorption  Profile  DPAL  Regime  (3/4  View) 


The  input  pump  was  matched  pretty  close  in  bandwidth  to  the  cross  section  at  800 
torr.  At  the  entrance  face  (0  cm)  of  the  cell  the  shape  of  the  input  pump  is  an  obvious 
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Gaussian.  As  the  power  is  absorbed  a  linear  decay  is  observed  at  every  frequency.  Right 
before  termination  however  the  profile  loses  its  linear  dependence.  This  suggest  that 
perhaps  a  linear  approximation  can  be  used  to  determine  a  cell  length  of  a  DPAL 
absorbing  system.  The  linear  nature  of  this  region  using  equation  1.89  is: 


7(x,v)  =  /0(v)  + 


2  Nxcj3  ,  u  y )  4,  u  y )  •  (  - 1  +  4  Kjl  h  ) 

(3  +  2k32T2A3]T3) 


(1.100) 


From  inspection  of  Figures  21  and  22  it’s  hard  to  discern  when  the  profile  ceases  to  be 
linear  and  how  fast  or  slow  it  transitions  from  being  linear  to  exponential.  Eventually  as 
intensity  is  absorbed  throughout  the  cell  at  some  point  that  intensity  will  be  at  or  below 
saturation  intensity  and  that  is  why  the  profile  rolls  off  as  it  nears  zero. 

To  see  this  more  closely  observe  the  profile  at  one  frequency,  in  this  case  the  4  to 
3  hyperfine  transition,  at  various  input  powers  of  Isat  plotted  against  the  linear 
approximation.  Figure  22  below  shows  such  a  plot.  The  four  sets  of  curves  shows  are, 
starting  from  highest  to  lowest,  input  powers  of  7,  5,  3,  and  1  times  of  Isat  labeled  a,  b,  c, 
and  d  respectively.  The  trend  is  clear,  the  higher  the  initial  pump  intensity  is  compared  to 
that  of  Isat,  the  more  the  absorption  profile  follows  a  linear  trend.  Even  at  seven  times  the 
saturation  intensity  the  profile  deviates  from  linear  behavior  rather  quickly.  At  about  'A 
centimeter  into  the  cell,  exponential  behavior  begins  to  be  mixed  into  with  the  linear 
behavior.  It  is  obvious  that  the  linear  approximation  fails  to  predict  at  what  distance  into 
the  cell  the  input  is  attenuated  to  zero,  by  more  than  50%,  when  the  pump  intensity  starts 
near  Isat. 
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2D  Absorption  Profile  of  DPAL  Regime  at  Various  Iq 
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Figure  23.  2D  Absorption  Profiles  of  DPAL  Regime  at  Various  Io. 

Differential  Transmittance  and  the  Bleached  Wave 

As  stated  already  another  quantity  can  be  explored  using  this  model  termed 
differential  transmittance.  This  dTransmittance  profile  is  produced  by  taking  the 
derivative  of  the  absorption  profile  with  respect  to  the  input  pump.  Such  a  profile  will  be 
able  to  tell  us  where,  if  added  to  the  system  after  absorbing  conditions  have  been 
established,  one  photon  will  leave  the  area  of  transparency  and  hit  a  saturation  curtain  at  a 
particular  frequency.  The  importance  of  such  profile  lies  in  the  fact  that  if  the  laser  is  to 
operate  at  peak  efficiency  one  must  make  sure  that  no  saturation  curtain  exists  in  the  cell. 
In  essence  this  area  of  transparency  is  called  a  bleached  wave  and  this  dTransmittance 
profile  is  a  visual  representation  of  this  wave. 

The  differential  transmittance  profile  for  the  Beer’s  Law  regime  is  rather 
uninteresting  and  is  shown  in  Figure  23  along  with  its  transmittance  plot  in  Figure  24. 


72 


Once  again  in  Figure  24  we  see  the  transmission  window  a  bandwidth  of  about  7000 
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Figure  24.  dTransmittance  Profile  Beer’s  Law  Regime 
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Figure  25.  Transmittance  Profile  Beer’s  Law  Regime 
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Figure  26.  dTransmittance  Profile  DPAL  Regime 
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Figure  27.  dTransmittance  Profile  DPAL  Regime 


The  dTransmittance  profiles  for  the  DPAL  regime  are  very  revealing  and  are  displayed  in 


Figures  25  and  26.  As  one  can  see,  after  absorbing  conditions  are  set,  if  one  photon  was 
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sent  through  the  entrance  of  the  cell  that  did  not  add  any  significant  intensity  to  the  pump 
it  would  travel  unobstructed  until  it  would  hit  an  opaque.  It  is  clear  that  this  profile  has  an 
abrupt  crash  where  the  transparency/opaque  interface  occurs,  hence  the  bleached  wave 
nomenclature.  In  the  case  for  these  DPAL  system  parameters,  in  order  to  maximize 
absorption  throughout  the  cell,  the  length  would  need  to  be  around  6  cm  in  length. 

This  being  an  inherently  numerical  model  the  question  arises  is  there  any  quick 
analytic  determination  to  the  length  of  the  bleached  wave?  Turning  to  the  linear 
approximation  Figures  27  and  28  examine  the  performance  of  this  approximation  to 
bleached  wave  characteristics. 
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2D  Absorption  Profile  of  DPAL  Regime  at  /o=7*Isat 
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Figure  28.  Linear  Approximation  Performance  at  a.)  7xlsat  and  b.)  5xlsat 
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2D  Absorption  Profile  of  DPAL  Regime  at  /o=3*Isat 
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Figure  29.  Linear  Approximation  Performance  at  a.)  3xlsat  and  b.)  Isat 
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The  same  single  frequency  analysis  is  done  again,  but  this  time  with  dTransmittance  and 


linear  approximation  curves  overlaid.  A  couple  of  bleached  wave  characteristics  can  be 
calculated  from  the  model:  distance  at  which  dT  equals  90%,  distance  at  which  it  equals 
10%,  and  the  distance  it  takes  for  the  wave  to  crash  from  90%  to  10%  normalized  to  the 
ideal  cell  length.  Also  the  value  at  which  the  linear  approximation  falls  to  zero  as  well  as 
what  is  the  dT  value  at  this  null  length  is  tabulated  in  Table  13. 


Table  14.  Linear  Approximation  and  dTransmittance  Performance  Numbers 


Pump 
Intensity 
(x  Isat) 

LinearZnull 

dT  (znull) 

90% 

(cm) 

10% 

(cm) 

90-10 

(cm) 

90-10 

(Normalized) 

1 

0.56 

68.32% 

0.24 

1.24 

1.00 

0.51 

3 

1.71 

69.58% 

1.04 

2.56 

1.52 

0.51 

5 

2.86 

71.17% 

2.02 

3.80 

1.78 

0.37 

7 

4.00 

72.36% 

3.07 

5.01 

1.94 

0.33 

As  discussed  before  the  linear  approximation  is  a  poor  representation  of  the  actual 
absorption  profile  at  input  pump  intensities  around  Isat  and  a  decent  representation  when 
intensities  are  way  above  Isat.  The  remarkable  aspect  of  the  linear  approximation, 
however,  lies  in  what  it  predicts  the  dTransmittance  is.  No  matter  what  the  input  intensity 
is  a  good  length  to  construct  an  absorbing  cell  is  the  distance  predicted  by  z  null.  It 
appears  that  about  70  of  100  single  photons  will  pass  through  without  being  absorbed. 
Another  aspect  that  can  be  discerned  is  how  hard  the  bleached  wave  crashes.  The 
distance  at  which  the  differential  transmittance  decreases  from  90%  to  10%  normalized  to 
its  cell  length  gives  an  indication  for  comparison.  The  lower  this  ratio  is  the  more  defined 
the  bleached  wave’s  front  is.  Visually  this  is  seen  in  Figures  27a  and  28b.  The  wave  in 
27a  falls  off  dramatically  where  as  the  wave  in  28b  is  more  of  a  gentle  slope. 
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Quasi  Two  Level  Analytic  Comparisons 

Quasi  Two  Level  Analytic  Solution 

An  analytic  solution  to  the  rate  equations  and  the  differential  law  of  absorption 
can  be  made  if  one  assumes  that  the  spin  orbit  relaxation  between  the  P3/2  and  P1/2  states 
is  the  fastest  process  compared  to  stimulation  and  emission.  Essentially  then  the  two 
states  cannot  be  distinguished  and  the  picture  of  the  system  becomes  that  of  a  quasi-two 
level  absorption  system.  The  two  upper  levels  are  statistically  related  by  Boltzmann 
statistics  as: 


/ 


1 


1  +  —  •  Exp(-AE32  /  kT) 
g 2 


(1.101) 


As  such  the  population  of  the  statistically  bound  upper  state  becomes  n2=  f  #  and 

n3  =  (1  -  f)nLI  and  the  ground  state  population  is  simply  labeled  n,  .  Combined  with  the 

same  constraint  that  the  total  population  does  not  increase  or  decrease  the  rate  equations 
for  the  quasi  two  level  system  become: 


dn‘'  ~  Ip  - g3r ((! ~ ZK  ~ 2nL ) +  ^/+,+ 1  ^  o~2i (/ 9  ~  ^)+— 


dt  hv. 


hvr 


dlp  _ 

dz  ±p  ~31 

^r=°2i w%-  5)  1 

dz 

=  ti)  L 

dz 


(1.102) 


79 


where  z  is  the  effective  lifetime  of  the  combined  upper  states  and  1+  and  L  are  the  right 
and  left  running  intracavity  laser  intensities.  A  similar  game  can  be  played  in  which  the 
rates  in  equation  1.104  are  set  to  steady  state  conditions  and  from  there  the  population 
inversion  is  calculated.  This  population  inversion  is  then  inserted  into  the  differential  law 
of  absorption  and  what  is  obtained  is  a  quasi  two  level  description  of  pump  absorption: 


1  dl, 


Ip  dz 


P  _ 


=  -cr3i(v> 


2+77(3/-!) 


v4(A> 


l  +  (3 -/)■ 


+  77(/  +  l)- 


(1.103) 


To  compare  with  the  rate  equation  analysis  of  this  model  the  lasing  intensity  in  equation 
1.105  is  set  to  zero  and  then  the  differential  equation  is  solved  for  simply  an  absorbing 
case.  Once  accomplished  the  same  analysis  done  for  the  rate  equation  analysis  (REA) 
method  is  done  for  the  quasi  two  level  absorption  analysis  (QTLA)  and  the  results 
compared.  Figures  29  and  30  display  both  the  absorption  profile  and  differential 
transmittance  profiles  in  space  and  in  frequency: 
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Figure  30.  QTLA  Profile 


dTransmiltance  (QTLA) 


-20000  MHz 


Figure  31.  QTLA  dTransmittance  Profile 


It  is  obvious  that  in  comparison  to  Figures  21  and  25  the  two  solutions  show  similar 
features.  Both  exhibit  a  linear  drop  off  in  intensity  and  an  exponential  fall  off  as  intensity 
nears  zero.  The  bleached  wave  pictures  look  similar  as  well.  Transparency  of  the  cell 
crashing  into  an  opaque  curtain,  however  for  the  same  parameters  it  appears  that  the 


bleached  wave  does  not  penetrate  as  far  into  the  cell  as  in  the  rate  equation  analysis.  To 
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further  compare  the  two  methods  Figures  31-34  displays  at  different  values  of  the 
saturation  intensity  the  single  frequency  profile  from  the  QTLA  analysis  versus  the  linear 
approximation  of  the  rate  equation  analysis.  Whereas  the  linear  approximation  undershot 
the  REA  curves  it  tends  to  be  right  on  top  of  the  QTLA  curves.  The  following  four 
figures  show  bleached  wave  characteristics  and  linear  approximation  performance  values. 
Table  14  is  augmented  with  these  new  performance  values  in  Table  15  for  ease  of 
comparison. 
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Figure  32.  Linear  Approximation  Performance  against  QTLA  at  7x  Isat 
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Absorption  Profile  QTLA  at  7o=5*Isat 
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Figure  33.  Linear  Approximation  Performance  against  QTLA  at  5x  Isat 
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Figure  34.  Linear  Approximation  Performance  against  QTLA  at  3x  Isat 
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Absorption  Profile  QTLA  at  /q=  1  *Isat 
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Figure  35.  Linear  Approximation  Performance  against  QTLA  at  lx  Isat 


Table  15.  REA  and  QTLA  Linear  Approx  dTransmittance  Numbers 


Pump 
Intensity 
(x  Isat) 

LinearZnull 

dT  (znull) 

90% 

(cm) 

10% 

(cm) 

90-10 

(cm) 

90-10 

(Normalized) 

1 

0.56 

68.32% 

0.24 

1.24 

1.00 

0.51 

3 

1.71 

69.58% 

1.04 

2.56 

1.52 

0.51 

5 

2.86 

71.17% 

2.02 

3.80 

1.78 

0.37 

7 

4.00 

72.36% 

3.07 

5.01 

1.94 

0.33 

1 

0.56 

63.50% 

0.21 

1.17 

0.97 

0.49 

3 

1.71 

59.54% 

0.91 

2.38 

1.47 

0.49 

5 

2.86 

56.26% 

1.78 

3.51 

1.73 

0.41 

7 

4.00 

52.24% 

2.72 

4.61 

1.89 

0.35 

Quasi  Two  Level  Absorption  vs.  Lasing 

Once  more  we  repeat  the  same  type  of  analysis  as  done  with  the  previous  two 
sections  using  equation  1.105  with  all  processes  active.  The  shape  and  function  of  the 
output  laser  intensity  is  highly  dependent  on  the  length  of  the  cell.  The  QTL  regime  uses 
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a  longitudinal  geometry  with  key  laser  parameters  being  a  function  of  cell  length  and  the 
reflectivity  of  the  output  coupler.  Also  for  lasing  to  take  place  a  certain  threshold  must  be 
crossed  or  else  the  laser  will  not  run.  The  following  quantities  have  been  derived 
following  the  QTL  approximation  and  derivation. 


Effective  Absorbance: 


/f  =  (3/zl) 

(/  +  1)  31  '0t  S 


(3~/) 

(/  + 1) 


where  /  is  the  cell  length  and  r  is  the  output  coupling  parameter. 


(1.104) 


Laser  Threshold  Intensity: 


4  = 


(1-^X1+/) 


31 

y^2i 


Ln 


+  cr3l  •  ntol  -lg 


(1.105) 


Slope  Efficiency: 


“I  slope 
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(1.106) 


Output  Laser  Intensity: 


Ioul=VsIope(  4(°)“4) 


(1.107) 
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In  order  to  setup  a  lasing  scenario  first  a  2  cm  cell  was  chosen  and  the  an  input 
power  was  picked,  70  W/cm2,  such  that  the  cell  was  nearly  completely  bleached.  The 
output  coupling  constant  was  chosen  to  be  0.5.  The  rest  of  the  system  parameters  were 
identical  to  the  DPAL  regime  parameters.  Based  on  the  system  parameters  the  following 
were  established  for  lasing  threshold  as  well  as  the  laser  output  power  in  frequency  space. 
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Figure  36. 1  Threshold 
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Figure  37.  Laser  Output 
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By  examining  equation  1.105  it’s  obvious  that  any  difference  between  the  absorption 
profile  with  and  without  lasing  will  occur  between  -4000  and  5800  MHz  from  inspecting 
Figure  35.  Equation  1.109  explains  that  above  threshold  the  output  intensity  varies 
linearly  with  respect  to  the  input  pump  intensity.  Taking  these  facts  into  consideration  we 
look  at  the  absorption  profiles  of  both  cases  at  zero  frequency  for  the  following  input 
intensities:  56  (threshold),  70,  140,  and  210  W/cm2. 
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Difference  of  Profiles  w/o  Lasing  and  w /  Lasing 
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Figure  39.  Difference  of  Profiles  w/o  Lasing  and  w/  Lasing  at  I0=70 


88 


Difference  of  Profiles  w/o  Lasing  and  w/  Lasing 
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Figure  40.  Difference  of  Profiles  w/o  Lasing  and  w/  Lasing  at  I0=140 
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Differeoce  of  Profiles  w/o  Lasing  and  w/  Lasing 


Figure  41.  Difference  of  Profiles  w/o  Lasing  and  w/  Lasing  at  I0=210 
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Figure  37  shows  no  marked  difference  in  absorption  profiles.  This  of  course  makes 
perfect  sense  because  the  input  intensity  is  at  the  required  threshold  for  lasing.  Figures 
38-40  however  have  the  necessary  intensity  to  achieve  lasing.  A  remarkable  result 
occurs:  no  matter  the  input  power  above  threshold  the  lasing  condition  seems  to  keep 
absorption  profile  fairly  constant.  The  top  curves  of  Figures  38-40  represent  the 
absorption  profile  without  lasing.  We  can  see  the  key  feature  of  linearity  when  the  input 
intensity  is  significantly  above  Isat  however  the  lasing  prevents  this  limit  from  being 
achieved.  This  implies  an  important  aspect  about  cell  size  and  lasing.  For  a  given  set  of 
parameters  the  calculated  intensity  for  lasing  threshold  must  at  minimum  be  able  to 
bleach  the  entire  cell  length.  If  this  is  achieved  Figure  38-40  imply  that  this  bleach  length 
will  remain  constant  no  matter  how  high  above  threshold  the  input  intensity  is. 
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IV.  Conclusions 


Simulation  Scenarios 

Beer’s  Law  Regime 

The  outputs  of  this  model  to  the  Beer’s  Law  regime  were  expected  since  the 
differential  law  of  absorption  has  an  analytical  solution  when  the  input  intensity  is  well 
below  saturation  intensity.  As  shown  in  Figure  20  however  with  the  rapid  absorption 
along  the  hyperfine  peaks  a  transmission  window  is  created  which  could  be  useful  for 
optical  systems.  The  advantages  of  such  a  system  is  first  the  bandwidth  of  the  window  is 
controlled  through  increasing  or  decreasing  the  system  pressure.  As  such  the  window 
becomes  a  user  controlled  variable.  Second  since  all  alkali  metals  have  similar  atomic 
structure,  mainly  the  doublet  absorption  cross  section,  different  center  frequencies  can  be 
chosen  by  using  a  different  alkali  metal.  As  such  one  could  engineer  a  variable  bandwidth 
transmission  window  whose  center  frequencies  are  those  of  each  alkali  metals  unique  fine 
structure  transition  frequency. 

It  is  also  not  too  surprising  that  the  dTransmittance  profile  looks  very  similar  to 
the  absorption  profile.  The  analytic  solution  to  Beer’s  Law  when  solved  and  the 
derivative  taken  with  respect  to  the  input  intensity  Io  simply  yields  another  exponential. 
The  advantage  of  the  Beer’s  Law  simplicity  lends  itself  to  be  an  easy  benchmark  to  test 
any  absorbing  model.  The  rate  equation  analysis  that  this  model  used  made  no  such 
Beer’s  Law  approximation,  however  it  is  expected  that  when  Beer’s  Law  conditions  hold 
the  full  solution  should  yield  the  analytic  Beer’s  Law  solution.  It  is  obvious  from  Figures 
20,  23,  and  24  that  the  rate  equation  analysis  yields  what  is  expected.  All  rate  equation 
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solutions  derived  or  modeled  should  always  be  first  tested  under  Beer’s  Law  conditions 
since  it  is  such  an  easy  benchmark  of  performance. 

DPAL  Regime 

The  DPAL  regime  took  a  look  at  the  other  end  of  the  spectrum:  high  system 
pressure  and  high  input  intensity.  No  such  doublet  absorption  was  observed  since  the 
hyperfine  peaks  were  pressure  broadened  into  one  peak.  By  far  pressure  has  the  most 
variable  effects  on  lineshape  broadening.  Natural  broadening  has  little  to  no  effect  in 
when  compared  to  Doppler  and  pressure  broadening.  Even  Doppler  broadening  has  its 
limitations  because  the  only  way  to  broaden  the  lineshape  is  to  increase  temperature  and 
if  the  temperature  is  continually  raised  eventually  one  might  start  melting  parts  of  their 
apparatus.  Not  to  mention  raising  system  temperature  for  the  purposes  of  broadening  a 
lineshape  could  introduce  thermal  regulation  issues  that  are  not  accounted  for  in  this 
model.  When  needed  pressure  should  always  be  the  chief  mechanism  to  broaden  any 
lineshape  in  such  an  absorbing  configuration. 

Upon  solving  the  differential  law  of  absorption  using  the  population  densities 
derived  from  the  rate  equations  while  the  input  power  was  above  Isat  proved  to  have  no 
analytical  solution.  In  fact  the  numerical  solution  involves  a  function  called  the  Product 
Log,  or  the  Lambert-W  function,  which  is  an  inherently  complex  function.  So  even  when 
modeling  this  solution  one  needs  to  understand  that  do  to  its  complex  nature  multiple 
solutions  exist.  Finding  the  right  solution  that  models  real  world  experiment,  if  they  exist, 
could  prove  to  be  a  significant  task  for  different  developments  of  rate  equations. 

In  the  upper  limit  that  the  input  intensity  was  significantly  higher  than  Isat  an 
approximation  to  the  differential  law  of  absorption  was  made  which  turned  out  to  be 
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analytic.  That  solution  was  the  linear  approximation.  How  high  the  input  has  to  be  above 
Isat  where  the  linear  approximation  becomes  an  accurate  account  of  the  absorption 
profile  simply  has  to  be  determined  experimentally.  From  Figure  22  the  actual  profile  is 
undoubtedly  a  mix  of  linear  and  exponential  character.  Even  at  an  input  intensity  of  7x 
Isat  the  exponential  quality  of  the  solution  takes  hold  rather  quickly.  Determining  where 
the  transition  from  linear  to  exponential  character  occurs  is  nearly  impossible  to  tell. 

Once  the  solution  to  the  saturation  law  using  the  rate  equation  analysis  was 
numerically  solved  taking  the  derivative  of  this  solution  with  respect  to  the  input  intensity 
led  to  a  new  quantity  termed  differential  transmittance.  The  differential  transmittance  is  a 
representation  for  the  bleached  wave.  The  wave  front  of  the  bleached  wave  marks  where 
the  region  of  transparency  ends  and  where  an  opaque  region  begins.  In  practice  an  actual 
system  would  need,  for  optimal  performance,  the  cell  length  to  match  the  length  of  the 
bleached  wave.  Any  cell  length  longer  than  this  wave  will  lead  to  a  hit  in  absorption 
performance  which  will  directly  degrade  the  laser  output  performance. 

The  null  distance  predicted  by  the  linear  approximation  while  poor  for  predicting 
where  the  input  is  attenuated  to  zero  gave  a  pretty  consistent  dTransmittance  prediction 
regardless  if  the  system  was  in  a  Beer’s  Law  or  DPAL  regime.  Table  14  shows  that  if  the 
linear  approximation  is  used  to  make  a  first  estimate  on  cell  length,  about  70%  of  single 
photons  will  bleach  through  the  cell.  Since  the  linear  approximation  has  such  a  simple 
form,  equation  1.102,  if  one  knows  beforehand  what  system  parameter  they  want  to  run, 
lifetimes,  and  kinetic  constants,  equation  1.102  simply  needs  to  be  inverted  and  the 
predicted  cell  length  will  be  a  very  good  estimate  to  run  an  absorbing  system.  Looking  at 
the  slope  of  equation  1.102, 
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(1.108) 


2-  N  z  lsat(v)  c7(u)  (l-  ^3^2-r2T3) 

(3  +  k32t2^A3lT3 

it  is  obvious  that  what  a  DPAL  system  is  most  sensitive  to  is  the  concentration  of  alkali. 
Barring  that  any  further  variations  due  to  pressure  and  temperature  will  be  contained  in 
how  much  lifetimes  and  the  saturation  intensity  vary.  Nonetheless  the  linear 
approximation  is  concluded  to  be  a  good  starting  point  in  creating  a  DPAL  absorbing 
cell. 

Comparison  to  Quasi  Two  Level  Analysis 

To  effectively  compare  the  rate  equation  model  with  the  quasi  two  level  model  we 
take  a  look  at  how  the  population  inversion  density  changes  with  initial  input  pressure  at 
particular  choices  of  system  pressure.  The  population  inversions  of  the  REA  model  and 
the  QTL  model  are 
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respectively.  When  there  is  no  pump  present  both  of  these  normalized  inversions  will 
evaluate  to  negative  two  and  as  the  pump  goes  to  infinity  these  values  will  evaluate  to 
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zero.  The  quasi  two  level  model  assume  beforehand  that  the  second  and  third  levels 
coupled.  This  statement  in  the  REA  model  amounts  to  a  very  high  system  pressure  for 
this  would  force  the  mixing  rate  to  be  the  fastest  rate  in  the  system. 

Figure  41  displays  a  set  of  these  curves  for  different  total  system  pressures.  At 
any  particular  input  power  we  see  that  it  is  harder  to  maintain  a  N3,  N1  population 
inversion.  For  instance  at  an  initial  pump  intensity  of  500  W/cnr  at  a  system  pressure  of 
400  torr  (curve  A)  the  value  of  the  inversion  is  close  to  zero  or  in  other  words  the 
population  in  the  first  and  third  states  are  comparable.  However  looking  at  curve  E, 
where  the  system  is  at  20000  torr,  we  find  more  population  in  the  first  state  than  the  third 
state.  The  reasoning  behind  this  is  that  at  high  system  pressures  the  mixing  rate  from  state 
3  to  state  2  is  faster  than  the  pumping  rate;  state  three  is  being  de-populated  and  it  is 
harder  to  maintain  the  N3,  N1  inversion. 

For  a  particular  pressure  when  the  normalized  population  inversions  are  compared 
between  the  two  models  we  expect  the  quasi  two  level  inversion  to  be  slightly  harder  to 
maintain  than  that  of  the  REA  model.  Figure  42  shows  the  two  models  compared  at 
system  pressures  of  400,  6000,  and  20000  torr.  In  each  case  the  QTF  inversion  maintains 
a  slight  edge  in  depopulating  its  third  level  to  the  second  level.  In  each  case  the  REA 
normalized  population  inversion  approaches  that  of  the  QTF,  but  never  surpasses  it.  This 
shows  that  as  system  pressure  is  increased  the  REA  model  will  predict  pump  attenuations 
closer  to  that  of  the  QTF  because  the  attenuation  depends  on  the  N3,  N1  population 
inversion.  The  REA  model  converges  to  the  QTF  model  under  QTF  assumptions. 
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Figure  42.  Normalized  Population  Inversions  at  various  system  pressures.  Curve  A: 
400  torr;  Curve  B:  2000  torr;  Curve  C:  6000;  Curve  D:  10000;  Curve  E:  20000 
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Figure  43.  Normalized  Population  Inversions  of  REA  model  vs  QTLA  model  at 
system  pressure  of  400,  6000,  and  20000  torr. 
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Comparison  of  QTL  Absorption  and  Lasing 

The  final  analysis  of  this  paper  was  to  compare  the  absorption  profiles  within  the 
QTL  analysis  between  lasing  off  and  lasing  on.  The  laser  parameters,  effective 
absorbance,  slope  efficiency,  and  threshold  were  all  solely  defined  around  system 
parameters  with  the  value  of  the  output  coupler  reflectivity  user  defined.  A  test  scenario 
was  the  set  up  in  which  a  2  cm  cell  was  chosen  and  the  system  parameters  the  same  as  the 
DPAL  regime.  Figures  36  display  such  a  comparison.  As  seen  in  Figure  37  the  initial 
input  intensity  was  at  threshold  so  no  lasing  effects  occurred  and  that  is  why  only  one 
profile  exists.  The  pictures  differ  drastically  in  Figures  38-40.  The  upper  curves,  which 
represent  no  lasing,  simply  blast  through  the  cell.  However  the  curves  when  lasing  is 
running  are  nearly  constant  regardless  of  what  the  input  power  is.  In  fact  they  look  really 
similar  to  the  profile  where  the  input  intensity  is  at  laser  threshold. 

To  make  sense  of  this  phenomenon  the  quasi  two  level  system  essentially  is 
describing  a  photon  engine.  A  photon  absorbs  energy  from  the  pump,  is  excited  to  the 
second  level,  and  then  returns  to  the  ground  state  by  stimulated  emission.  This  process 
happens  over  and  over  again  hence  the  photon  engine  moniker.  With  this  repeated 
process  pump  intensity  must  be  devoted  to  those  photons  that  have  run  this  cycle 
repeatedly.  This  is  why  the  absorption  profiles  with  lasing  on  in  Figures  38-40  remain 
nearly  the  same.  The  reason  why  they  are  not  the  same  as  the  ideal  profile  in  figure  37  is 
because  quantum  and  slope  efficiencies  will  never  be  100%.  The  further  the  laser 
quantum  efficiency  deviates  from  100%  the  less  likely  the  absorption  profile  will  mimic 
the  ideal  profile  in  which  the  input  intensity  is  threshold  intensity.  This  discovery  leads  to 
another  engineering  benchmark.  The  laser  threshold  value,  computed  from  system 
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parameters,  must  be  able  to  entirely  bleach  the  selected  cell  length.  If  this  is  not  the  case 
then  adequate  absorption  cannot  take  place  and  the  laser  will  not  run. 

Combining  this  result  with  the  linear  approximation  for  a  quasi  two  level  system 
we  have  a  benchmark: 

m  °  4,,^,)- ' 2  am) 

Once  a  desired  cell  length  is  chosen,  zq,  and  inserted  into  equation  1.113,  if  the  resultant 
answer  is  less  than  or  equal  to  zero  and  the  threshold  is  much  higher  than  the  saturation 
frequency,  so  that  the  linear  approximation  is  fairly  accurate,  the  laser  will  not  run. 

Recommendations  for  Future  Modeling  Work 

The  next  step  for  this  model  to  take  would  be  to  implement  lasing  from  a  Rate 
Equation  Analysis  standpoint.  This  would  involve  deriving  additional  rate  equations  that 
account  for  the  photons  being  created  along  the  lasing  D1  manifold  as  the  spin  orbit 
processes  regulate  population  from  the  D2  to  D1  manifold.  Solving  the  equations  would 
be  a  mere  exercise  in  Mathematica  modeling  skill,  but  inherent  in  the  lasing  rate 
equations  will  be  the  next  big  problem  for  this  model  to  solve:  what  are  the  ideal 
geometrical  configurations  that  would  maximize  laser  output?  One  could  vary  the  cell 
from  a  longitudinally  pumped  cell  to  a  transversely  pumped  cell.  Or  perhaps  the  cell 
could  be  longitudinally  pumped  but  lase  transversely?  The  size  and  shape  of  the  cavity 
would  force  different  cavity  modes  to  be  explored  and  analyzed  and  the  configurations 
are  practically  limitless.  So  by  establishing  two  to  three  different  geometrical 
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configurations  and  then  solving  pump  rates  and  lasing  rates  a  stronger  engineering  hold 
of  what  a  DPAL  system  should  look  like  will  come  more  to  light. 

Recommendations  for  Future  Experimental  Work 

The  robustness  of  this  model  must  be  tested  against  a  wide  range  of  scenarios. 
Two  particular  experiments  can  test  this  model  right  away.  The  first  experiment  would  be 
absorption  along  a  longitudinal  axis  of  an  absorbing  cell.  Multiple  test  scenarios  should 
be  run  at  a  wide  variety  of  temperatures,  pressures,  input  powers,  and  varying 
compositions  of  foreign  species  especially  those  scenarios  that  lead  to  a  bleached  wave 
being  formed.  This  way  profiles  can  be  compared  readily  to  those  outputted  by  this 
model. 

A  second  experiment  can  be  devised  to  test  the  bleached  wave  and  linear 
approximation.  First,  bleaching  conditions  should  be  established  based  of  a  cell  length 
predicted  by  the  linear  approximation.  Second,  once  bleaching  conditions  are  established 
a  secondary  pump  slight  detuned  from  the  primary  pump  can  be  used  at  low  power  to 
inject  photons  into  the  system  that  would  not  affect  the  primary  input  pump  power.  If 
these  low  power  photons  can  be  detected  on  the  output  end  of  the  cell  an  experimental 
picture  of  the  bleached  wave  can  be  detected  and  then  compared  to  the  theoretical 
dTransmission  value  predicted  from  the  linear  approximation.  Then  system  parameters 
can  be  slightly  adjusted  such  that  the  bleached  wave  creeps  back  into  the  cell  or 
completely  blasts  through  the  length.  The  same  low  power  photons  can  be  injected  to  see 
if  the  distance  of  the  bleached  wave  changes  as  predicted  by  this  model. 
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Appendix  A 


Development  of  Population  Inversion 
Rate  Equations  for  a  three  level  system: 

(2.1) 

(2.2) 

(2.3) 


dN, 


3  _ 


dt 


dN2 

dt 


N* 


r-IP(v)a3]  1  A 

y  /?Vp  T3 


+  N2k23  +  N3 


r  i p(y)o 3i  &  , A 

‘  *13 


V 


hvP  gl 


N3k22  +  N2 


-iL(y)< 7n  1 

hv,  T. 


+  A, 


g2 

hVL  g 


+  k* 


12 


dN  ( dN. '  ) 


dt 


-JL  + - 1 

^  dt  dt  j 


Substitute  the  total  number  of  alkali  atoms,  iV2  =  A'  -  A,  -  A), ,  into  equation  2.1. 
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Set  all  time  derivatives  to  zero  for  steady  state  conditions  and  rearrange  equations  2.2, 
2.3,  and  2.4  into  a  matrix  equation: 
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Multiply  both  sides  by  4>i^3i  'n  order  to  get  equation  2.3  in  terms  of  saturation 
intensities.  Use  Kramer’s  Rule  for  a  3x3  matrix  to  solve  for  N3,  N2,  and  Ni. 
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The  solution  to  each  component  becomes: 


101 


_  Det[S 3] 
~~  Det[S ] 

_ 

~~  De/[5] 
_  Det[Si  ] 
~~  De/tS'] 


(2.9) 


O' 

The  population  inversion,  N 3 - -TV ,  becomes,  with  lasing  set  to  IL  =  0,  equation  2.10 

Si 


below: 


2  +  ku t3  +  k2ikl2T2T3  +  2ki2k2iT2T2 ) 


(1  +  k^T2  +  (A:]3  +  k22(k}2  —  k22yi2  +  k^2k2 2t.,)t3)  H - (3  +  k^T^  +  2k21T1')A2{C2 

1 SP 


(2.10) 


102 


Bibliography 


Demtroder,  Wolfgang.  Laser  Spectroscopy:  Basic  Concepts  and  Instrumentation.  Berlin: 
Springer-Verlag,  2003. 

Gross,  Kevin  C.  Class  notes,  PHYS  661,  Atomic  and  Molecular  Spectroscopy. 
Department  of  Engineering  Physics,  Air  Force  Institute  of  Technology,  Wright- 
Paterson  AFB  OH,  Winter  Quarter  2008. 

Hager,  Gordon  D.  and  others.  “A  Quasi-Two  Fevel  Analytic  Model  for  End  Pumped 
Alkali  Metal  Vapor  Fasers”.  Department  of  Engineering  Physics,  Air  Force  Institute 
of  Technology,  Wright-Patterson  Air  Force  Base,  OH.  (Draft  May  2008). 

Hilborn,  Robert  C.  “Einstein  coefficients,  cross  section,  f  values,  dipole  moments,  and  all 
that,”  Am.  J.  Phys.,  50:  982-986  (1982).  (Revised  February,  2002). 

Konefal  Z.,  “Observation  of  collision  induced  processes  in  rubidium-ethane  vapour,” 
Optics  Communication,  164:95-105  (1999). 

2  2 

Krause,  F.,  “Collisional  Excitation  Transger  Between  the  P1/2  and  “P3/2  Fevels  in  Alkali 
Atoms,”  Applied  Optics,  Vol.  5,  No.  9:  1375-1382  (Sep  1966). 

Shankar,  Ramamurti.  Principles  of  Quantum  Mechanics  Second  Edition.  New  York: 
Plenum  Press,  1994. 

Steck,  Daniel  A,  “Cesium  D  Fine  Data,”  available  online  at  http://steck.us.alkalidata 
(revision  2.1,  September  2008). 

Steck,  Daniel  A,  “Rubidium  85  D  Fine  Data,”  available  online  at  http://steck.us.alkalidata 
(revision  0.2,  1  September  2008). 

Steck,  Daniel  A,  “Rubidium  87  D  Fine  Data,”  available  online  at  http://steck.us.alkalidata 
(revision  2.1,  1  September  2008). 

Verdyne,  Joseph  T.  Laser  Electronics  Second  Edition.  Englewood  Cliffs,  New  Jersey: 
Prentics-Hall,  1981 


103 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  074-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data 
sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other 
aspect  of  the  collection  of  information,  including  suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information 
Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other 
provision  of  law,  no  person  shall  be  subject  to  an  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 
PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


2.  REPORT  TYPE 


1.  REPORT  DATE  (DD-MM-YYYY) 

27-03-2009 


4.  TITLE  AND  SUBTITLE 


Master’s  Thesis 


3.  DATES  COVERED  (From  -  To) 

March  2008  -  March  2009 


5a.  CONTRACT  NUMBER 


Theoretical  Model  Analysis  of  Absorption  of  a  Three  level  5b-  grant  number 
Diode  Pumped  Alkali  Laser  5c.  program  element  number 


6.  AUTHOR(S) 

Lewis  II,  Charlton  D.,  1st  LT,  USAF 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAMES(S)  AND  ADDRESS(S) 

Air  Force  Institute  of  Technology 

Graduate  School  of  Engineering  and  Management  (AFIT/EN) 
2950  Hobson  Way,  Building  640 
WPAFB  OH  45433-8865 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

High  Energy  Laser  Joint  Task  Office 
901  University  Blvd  SE,  Suite  100 
Albuquerque,  NM  87106 
Ackerman,  Harro  (505)248-8108 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


13.  SUPPLEMENTARY  NOTES 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFIT/GAP/ENP/09-M07 


10.  SPONSOR/MONITOR’S 
ACRONYM(S) 

HEL-JTO 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


14.  ABSTRACT 

This  paper  models  the  absorption  phenomena  of  light  in  a  three  level  diode  pumped  alkali  laser 
system.  Specifically  this  model  calculates  for  a  user  defined  set  of  system  parameters  the 
attenuation  of  the  input  pump  beam  and  characteristics  of  the  bleached  wave.  Using  Wolfram’s 
Mathematica  6.0  software  all  necessary  physics  for  an  accurate  description  of  absorption  was 
modeled  from  first  principles:  energy  levels,  cross  sections,  spin-orbit  kinetic  processes, 
saturation  frequencies,  pump  attenuation,  and  differential  transmittance,  which  is  a  representation 
of  the  bleached  wave.  A  specific  DPAL  scenario  was  simulated,  455K  system  temperature,  alkali 
concentration  of  6.1  •  10  ,  and  a  system  pressure  of  200  torr  of  He  and  600  torr  of  Ethane.  For  a 
range  of  initial  input  intensities  the  linear  approximation  to  the  beam  attenuation  predicted 
consistently  a  differential  transmittance  value  of  70%.  It  was  concluded  that  the  linear 
approximation  is  a  good  indicator  of  the  distance  a  bleached  wave  penetrates  an  absorbing  cell. 
This  model  was  also  benched  marked  against  the  quasi  two  level  model.  In  the  limit  of  high 
system  pressure  the  simulated  model  converged  to  the  same  population  inversion  as  that  of  the 
quasi  two  level  regime.  Finally,  within  the  quasi  two  level  regime  a  closed  analytic  equation  was 
developed  to  describe  under  what  conditions  the  system  would  lase. 


15.  SUBJECT  TERMS 

DPAL,  Bleached  Wave,  Gas/Electric  Laser 


17.  LIMITATION 

18. 

OF 

NUMBER 

ABSTRACT 

OF 

PAGES 

UU 

115 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Glen  P.  Perram,  PhD 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

(937)  785-3636  x4504 

(glen.perram@afit.edu) 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 


